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We develop a systematic theory of microwave-induced oscillations in the magnetoresistivity of 
a two-dimensional electron gas, focusing on the regime of strongly overlapping Landau levels. At 
linear order in microwave power, two novel mechanisms of the oscillations ( "quadrupole" and "pho- 
tovoltaic") are identified, in addition to those studied before ("displacement" and "inelastic"). The 
quadrupole and photovoltaic mechanisms are shown to be the only ones that give rise to oscillations 
in the nondiagonal part of the photoconductivity tensor. In the diagonal part, the inelastic con- 
tribution dominates at moderate microwave power, while at elevated power the other mechanisms 
become relevant. We demonstrate the crucial role of feedback effects, which lead to a strong inter- 
play of the four mechanisms in the nonlinear photoresponse and yield, in particular, a nonmonotonic 
power dependence of the photoconductivity, narrowing of the magnetoresonances, and a nontrivial 
structure of the Hall photoresponse. At ultrahigh power, all effects related to the Landau quanti- 
zation decay due to a combination of the feedback and multiphoton effects, restoring the classical 
Drude conductivity. 

PACS numbers: 73.40.-c, 78.67.-n, 73.43.-f, 76.40.+b 



I. INTRODUCTION 

In recent years, the nonequilibrium properties of quan- 
tum Hall systems in a moderate perpendicular magnetic 
field B have become a subject of intense experimental 1-30 
and theoretical 31-68 research (for review see Refs. 69-72). 
Much attention has been attracted to the experimental 
discovery of a novel type of l/£?-periodic resistivity os- 
cillations which arise in these systems under microwave 
illumination. 1,2 Remarkably, it was demonstrated soon 
after the first experiments that at higher radiation power 
the minima of the oscillations evolve into novel "zero- 
resistance states" (ZRS). ,5 

As was shown in Ref. 31, the ZRS can be understood as 
a direct consequence of the oscillatory photoconductivity 
(OPC), provided the latter may become negative. In- 
dependently of the microscopic origin of the oscillations, 
the negative OPC signifies an instability leading to the 
formation of spontaneous-current domains, which in turn 
yields the vanishing of the observable resistance. 

Unlike the explanation of the ZRS, so far there 
has been no common agreement as to the micro- 
scopic origin of the OPC. The most frequently studied 
mechanism 32,33 ' 36 - 39 ' 44,45,47,48 ' 55 - 57 ' 61 - 65 ' 67 of the OPC, 
suggested long ago in Ref. 32, implies that the oscilla- 
tions occur due to complex scattering processes in which 
electrons simultaneously are scattered off impurities, ab- 
sorb (emit) microwave quanta, and are displaced along 
the applied dc electric field. We term this mechanism 
a "displacement" mechanism. A different mechanism, 
called here the "inelastic" mechanism, was recently pro- 
posed in Ref. 34 and studied in detail in Refs. 40,50 (sim- 
ilar ideas were also discussed in Ref. 7). In the inelastic 
mechanism, the magnetooscillations of the dc current are 
generated by a microwave- induced change of the isotropic 



time- independent part -Fbo (s) of the electron distribution 
function, 



f(e, <p, t) = F vn (e) exp(if(p + inuit) 



(1.1) 



where e is the electron energy, <p the angle of the quasi- 
classical momentum, u) the microwave frequency. 

The above mechanisms have much in common. In both 
of them, the OPC originates from oscillations of the den- 
sity of states (DOS) v{e) of disorder-broadened Landau 
levels with changing e. Due to the oscillatory i>(e), op- 
tical transitions in the system lead to a nonequilibrium 
correction to the distribution function / which oscillates 
with both e/uj c and lu/lo c , where lu c is the cyclotron fre- 
quency. In turn, the u> / w c -oscillations of / translate into 
the observed w/w c -oscillations of the magnetoresistivity. 

A crucial difference between the two mechanisms is the 
following. In the displacement mechanism, the radiation 
directly affects the first angular harmonic F\q(e) of the 
distribution function (and thus the dc current), whereas 
all effects of the ac and dc fields on the even angular har- 
monics are neglected. By contrast, in the inelastic mech- 
anism the ac field influences the isotropic part Fqq(e) 
only. A calculation of the dc response in the resulting 
state with a modified Fqq(e) yields the OPC. While in 
the first case the effect on the current is gained directly in 
every microwave-assisted impurity-induced electron scat- 
tering event, in the latter case the effect is accumulated 
over a long period of time during which the electron dif- 
fuses in the field of impurities until it experiences inelas- 
tic scattering. At experimentally relevant temperatures 
T~ IK, the inelastic relaxation is governed by electron- 
electron scatterings 40,50 with the inelastic relaxation time 
r in cxT- 2 . 

In the linear regime with respect to the microwave 
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power, the contributions to the OPC generated by the 
displacement and inelastic mechanisms sum up indepen- 
dently. A comparison shows [see Eqs. (16) and (17) of 
Rcf. 50 and Eqs. (6) and (11) of Rcf . 39] that both mech- 
anisms produce the OPC with the experimentally ob- 
served phase, period, and S-damping of the oscillations. 
However, the slow inelastic relaxation rate (compared to 
the rate of electron collisions with impurities, r -1 ) makes 
the inelastic contribution larger by a factor Ti n /r q 3> 1. 
According to the calculation in Ref. 50, Ti n /r q <~ 100 un- 
der the experimental conditions. 73 Apart from the mag- 
nitude of the effect, the two contributions are qualita- 
tively different in their dependence on T and polariza- 
tion of the radiation. In accord with the experiments, 
the inelastic contribution decreases as t;„ oc T~ 2 with in- 
creasing T and does not depend on the direction of linear 
polarization of the microwave field. By contrast, the dis- 
placement mechanism yields a T independent contribu- 
tion which depends essentially on the relative orientation 
of the microwave and dc fields. Despite the above argu- 
ments, in a number of recent WO rks 44,45 ' 47,48 ' 55 ~ 57 ' 61 ~ 65 ' 67 
on this subject the inelastic mechanism has not been 
taken into account. 

In this work, we perform a systematic study of all rele- 
vant contributions to the OPC both in the linear regime 
and in several nonlinear regimes which emerge with in- 
creasing microwave power. On the experimental side, this 
work has been motivated by a number of interesting non- 
linear phenomena observed in recent experiments. 25-30 
From the theoretical point of view, a unified approach 
to the problem is necessary since at high microwave 
power the interplay of the inelastic and displacement 
mechanisms becomes essentially important. We show 
that the different mechanisms of the OPC strongly af- 
fect each other due to feedback effects already at mod- 
erate power levels which are easily accessible in the ex- 
periment. In particular, while in the linear regime the 
displacement and inelastic contributions to the OPC add 
up in phase, in the regime of a "saturated inelastic contri- 
bution" (SIC) the displacement contribution changes its 
sign, which leads to a nonmonotonic behavior of the pho- 
toresponse (still dominated by the inelastic contribution) 
with growing microwave power. In the ultrahigh power 
regime, all contributions to the dc conductivity that are 
related to the Landau quantization decrease with increas- 
ing power due to a competition between the feedback and 
multiphoton effects. As a result, the oscillations of the 
photoresponse vanish and the classical Drude conductiv- 
ity tensor is restored in the limit of high radiation power. 

On top of the nonlinear interplay between the inelastic 
and displacement mechanisms, we find two novel mech- 
anisms leading to the OPC, "quadrupole" and "photo- 
voltaic". In the quadrupole mechanism, the microwave 
radiation leads to excitation of the second angular har- 
monic F 2 o of the distribution function. The dc response 
in the state with nonzero F 20 yields an oscillatory contri- 
bution to the Hall part of the photoconductivity tensor 
which violates Onsager symmetry. In the photovoltaic 



mechanism, a combined action of the microwave and 
dc fields produces non-zero temporal harmonics F21 and 
-F01, so that the distribution function (1.1) acquires os- 
cillatory time dependence. The ac response in the result- 
ing state with excited F 2 i and F i contributes to both 
the longitudinal and Hall parts of the OPC. Provided 
7m/fq S> 1, the inelastic mechanism still gives the dom- 
inant contribution to the diagonal part of the photocon- 
ductivity tensor. However, the quadrupole and photo- 
voltaic mechanisms are the only ones yielding oscillatory 
corrections to the Hall part. 

The paper is organized as follows. In Sec. II we for- 
mulate the problem and overview the main steps in the 
derivation of the quantum kinetic equation obtained in 
Ref. 39. Using this kinetic equation, in Sec. Ill we put 
forward a general classification of contributions to the 
photoconductivity at first order in the microwave power. 
We discuss here four distinctly different mechanisms of 
the OPC. In Sec. IV we write down the kinetic equation 
for the case of strongly overlapping Landau levels, which 
we focus on hereafter. In Sec. V we calculate the lin- 
ear photoconductivity. In Sees. VI and VII we analyze 
several nonlinear regimes the system passes through with 
increasing microwave power. Namely, in Sec. VI we study 
the SIC regime at moderate power, driven by feedback 
effects. In Sec. VII we turn to the interplay of the feed- 
back and multiphoton effects in the regime of ultrahigh 
microwave power. In Sec. VIII we consider the possibility 
of experimental observation of the nonlinear effects in the 
magnetoresistivity measurements. In Sec. IX the range of 
applicability of the theory is discussed. The main results 
are summarized in Sec. X. 



II. 2D ELECTRON GAS IN A CLASSICALLY 
STRONG MAGNETIC FIELD UNDER 
MICROWAVE IRRADIATION 

Our theory is based on the approach to the problem 
of kinetics of a 2D electron gas (2DEG) in the presence 
of magnetic field and radiation that was developed by 
Vavilov and Aleiner. 39 This approach was used for a sys- 
tematic study of the displacement mechanism of the OPC 
in Ref. 39. It was further developed to treat the inelastic 
mechanism in Ref. 50. Based on controllable approxi- 
mations, the approach allows us to analytically describe 
the nonlinear behavior of the system under intense mi- 
crowave radiation. 



A. Model and parameters 

Let us first specify the model. The relation between 
the main parameters of our theory, which is also satisfied 
for the characteristic parameters of the 2DEG studied in 
the experiments on the OPC, is as follows: 

e F >T,w,w c ,r q - 1 >r t -; 1 ,T i ; 1 . (2.1) 



3 



Here sp is the Fermi energy, T the temperature, u> the 
microwave frequency, ui c the cyclotron frequency, r q and 
Ttr the quantum (single-particle) and transport disorder- 
induced scattering times, respectively, and r; n the inelas- 
tic relaxation time. The conditions (2.1) imply: 

• quasiclassical kinetics: the 2D electron states that 
contribute to transport belong to high Landau lev- 
els in a narrow energy strip of width T around the 
Fermi level, ep 3> T,lu c ; 

• predominantly small-angle scattering (smooth dis- 
order), T tr > T q ; 

• classically strong magnetic field, w c Tt r 3> 1; 

• pronounced oscillations in the DOS due to the Lan- 



dau quantization, uj 



c'q 



1. 



In the high-mobility structures used in the experiment 
a smooth random potential in the plane of the 2DEG 
was created by remote donors separated from the 2DEG 
by a spacer of width £/2 3> kp , where kp is the Fermi 
momentum. The Fourier transform of the correlation 
function of the random potential, 



W(q) = W(0)e 



(2.2) 



falls off exponentially for momentum transfers q 3> 
leading to T q /r tr = {kp^)^ 2 <C 1, where the quantum 
and transport scattering times in zero magnetic field are 
given by 



dq W(q) 
2tt 2v F 



OC 

Tq J 



J_ = Jdq W(g) f g \ 2 
Ttr J 2tt 2vp \2kp) 



(2.3) 
(2.4) 



Effects of the smooth disorder on the spectral and trans- 
port properties of 2D electrons occupying high Landau 
levels can be properly described by the self-consistent 
Born approximation (SCBA), 74 provided that £ does not 
exceed the magnetic length 75 (which is the case for the 
relevant magnetic fields in most of the experiments). 



B. Quantum Boltzmann equation 

Using the SCBA under the conditions listed in 
Eq. (2.1), Vavilov and Alciner derived the quantum 
Boltzmann equation (QBE), 39 which is the starting point 
of our calculation. The derivation of the QBE includes a 
few key steps which we highlight below. 

• Moving coordinate frame. 

The initial problem of electron kinetics in the quan- 
tizing magnetic field and the external (microwave 



+ dc) electric field E(t) in the presence of static im- 
purities is reduced to a problem of kinetics in the 
presence of "dynamic" impurities, whose potential 
is t dependent, by changing to a moving coordinate 
frame r — > r — C(i), where £(f) obeys 



dtC(t) 



d. 



UJ C £ 



eE(i) 



m e 



= 1 



(2.5) 



All effects of the external electric field are now 
included in the time dependence of £(t). The 
transformation (2.5), which is particularly con- 
venient for treating the external electric field 
nonperturbatively, is equivalent to the transfor- 
mation to Floquet states used in a number of 
works 44 - 45 ' 47 - 55 ' 57 ' 62 ' 63 - 67 on the OPC. 

• Keldysh equations within the SCBA. 

Within the SCBA, the Green's functions G {a) and 
the self-energies [(a) — > R,A,K refers to the 
retarded, advanced, and Keldysh components, re- 
spectively] are related to each other as 



•<<*) 
J 21 



w? w(q)e 

e iqr q{o) e -iqr 



21 



(2.6) 



where the subscript (21) denotes the times t?, and 
t\ on the Keldysh contour and £ 21 = C(^) — C(^i)- 

Quasiclassical approximation. 

Most importantly, using the conditions (2.1), Vav- 
ilov and Alciner reduced Eq. (2.6) to a simpler qua- 
siclassical equation in the "act ion- angle" represen- 
tation: 



(«)/ 



9& \^) 



lU) r 



k 



G { 2 ?(h + k; $) , 



(2.7) 
(2.8) 



where the operators h and ip are canonically con- 
jugated, [n, if] = —i. The eigenvalues of h and tp 
are the Landau level index and the angle coordi- 
nate of the momentum, respectively. The explicit 
expression for the linear operator K is given below 
in Sec. IIC. 

A significant difference between Eqs. (2.6) and (2.7) is 
that for high Landau levels the self-energy becomes h 
independent. It follows that the distribution function /, 
defined in the usual way by 



G R -G A -G K = 2 [ G R f - fG A } , 



(2.9) 



commutes with (p. Accordingly, the operator <p, which 
enters Eq. (2.7) and the impurity collision integral 

i St im {/} - t R f - f± A + \{t K + t A - t R ) , (2.10) 
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can be treated as a c-rmmber. Substitution of Eqs. (2.7) 
and (2.8) into Eqs. (2.9) and (2.10) leads to the following 
quantum kinetic equation: 



£{/} 

£{/}21 
St; m {/}21 



St im {/} (2.11) 
(9t+w c ^)/2i-St in {/}2i , (2.12) 

£2l(023/31 - Mg&) 



-h 



Here 



/3l£23523 + hsfcsigtl 



t=(h+t 2 )/2 



(2.13) 



(2.14) 



is the "center-of-mass" time. The inelastic collision in- 
tegral St; n {/}21 accounts for electron-electron scattering 
and for coupling to a thermal (phonon) bath. 

The spectrum of the problem is found from the coupled 
set of Eq. (2.7) and the equation of motion 



d \ - 

i— - iv c hj G^(n, 0) 



5{t 2 -h) 

2tt 



(2.15) 



Altogether, Eqs. (2.7), (2.11), and (2.15) determine both 
the spectrum and the distribution function. The dc cur- 
rent in the presence of microwave radiation is then given 
according to Eq. (3.30) of Ref. 39 by 



j — 3d + ind > 



3nd = en e dtC(t) , 

/ COS If 



(2.16) 
(2.17) 



2Re dr gK t+T f t+T!l 



(2.18) 



where j„d and jd stand for the nondissipative and dis- 
sipative components of the current, respectively, and n e 
is the electron concentration. Here and below, the bar 
denotes time averaging over the period of the microwave 
field in the stationary state, the angular brackets denote 
averaging over the angle <p. 



C. Kernel of the QBE 

To complete the formulation of the model, we now 
specify the kernel JC in Eqs. (2.7) and (2.11). Accord- 
ing to Ref. 39, it is can be represented as 



with all effects of the external electric field E(£) being 
incorporated in the function X 2 i(ip) <x C21 ^ E. Specifi- 
cally, we parametrize E(f) as 



E(t) = E dc 



V2 



Rc 



IS- 



+ 5+ 
— is + 



-iut 



(2.20) 



where Ed c is the dc electric field, the amplitude of the 
microwave field E w , and real s± with s+ + s 2 _ = 1 define 
elliptical polarization of E w with the main axes along x 
and y. In particular, s + — 1 (s_ = 1) corresponds to 
passive (active) circular polarization, s + = s_ = l/y/2 
to linear polarization along the x-direction (<p = 0). 
Using Eqs. (2.5) and (2.20), X = X 2 i(ip) is written as 



X 



C(*2)-C(*i) f-sinp 

£ \ COS if 

Xdc + Xuj , 



where 



<%dc = 
Xu; 



■ ( COS If 

— Edc 

Tt r UJ C \ Sin (fi 



U)t- \ -\ 

— £± cos(ip ± Ujt) , 

1 ± 

t- = t 2 — h, t is given by Eq. (2.14), and 
£± = s± 



sin ■ 



/2r q eE u v F 

T tT LU(UJ C ±LU) 



(2.21) 

(2.22) 
(2.23) 

(2.24) 



After a series expansion in 1/JcfS, "C 1 the kernel (2.19) 
acquires the final form we use below throughout the pa- 
per: 



K 

r q /Cj_ 



Tq/Cj 



T JC 



r" 1 + JC ± + JCj + JC V + o[(k F zy 
X 2 



X 




{l + X^f 



Z 9 * (1 + **)3 ^ 



(2.25) 
(2.26) 



(2.27) 
(2.28) 



Equations (2.25)^(2.28) correspond to Eqs. (3.48) and 
(3.49) of Ref. 39. We now turn to a systematic analysis of 
the QBE, starting with a general classification of various 
contributions to the OPC at order 0{E^). 



III. MECHANISMS OF THE OSCILLATORY 
PHOTOCONDUCTIVITY 



^ = ^/ dxexp G§e) 



x e-l x l- ix * 21 ^exp 



(2.19) 



In this section we present a general solution to the 
QBE (2.11) to first order in both the microwave power 
and the dc electric field. The case of strongly overlapped 
Landau levels, u-> c T q < 1, on which we focus in this pa- 
per will be considered at this order in more detail in 
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Sec. V. Throughout the paper we neglect any effect of 
the microwave radiation on the functions g R ( A ^ which 
determine the DOS: 

v(s,t) = 2u Q Re J^J dt,e lst -g^ t _^ t _ t _ /2 ^) . 

(3-1) 

This approximation is well justified in the limit ut c T q <C 1. 
However, at uj c r q <~ 1 the microwave and dc fields 
may lead to a pronounced modification of the DOS (see 
Sees. V and VI of Ref. 39). A manifestation of the effect 
of the external electric fields on the DOS in the photore- 
sponse will be discussed elsewhere. 78 

A. Classical limit 

To better understand the role of different terms in the 
kernel (2.25), it is instructive to look first at the QBE 
in the classical (with respect to the magnetic field) limit 
w c r q — > 0. In this limit, the DOS v = v a = m/2ir is 
constant, the functions g^\ = —g 2 \ = <K*2 — *i)/2, and 
the collision integral St; m {/}2i = £21/21 — /21-K0, where 
Kq = t^ 1 results from the action of the operator ft tt 
on unity. After the Wigner transformation, Eq. (2.11) 
reduces to 

(d t + iu c d v )f~st in {f} 

= ~ + J §K(n).f(e-n). (3.2) 

In the absence of electric fields (X = 0), the r.h.s. of 
Eq. (3.2) is written as K v f = fl!/, which describes 
a diffusive relaxation of the electron momentum in a 
smooth random potential. This part of the kernel, which 
is the smallest one in the parameter l/fcp£ <C 1, plays no 
role in this paper. 

The linear response to a small electric field X <C 1 is 
fully governed by the part Kj of the kernel, which is odd 
in E(i). Substituting JCj into Eq. (3.2) immediately gives 
the Drude formula for the conductivity. The term /Cj_ , of 
zero order in <C 1, is even in the electric field. It 

modifies the isotropic part of the distribution function, 
which leads to the effect of heating. Also, JC± is respon- 
sible for the electric-field-induced changes in higher even 
angular harmonics of / (as we show below, this yields 
additional contributions to the photoconductivity). 

It is worth noting that under the standard assumption 
of an energy-independent Tt r in the degenerate Fermi gas, 
all nonlinear (in Ed c and E w ) corrections to the Drude 
formula vanish to zero at the classical level. 39 ' 41 In or- 
der to get a finite dc photoresponse, one should then 
include an energy variation of r tr around the Fermi level, 
which is, however, still insufficient to get the oscillatory 
photoresponse. 41 As far as the OPC is concerned, at the 
classical level it can only come from the non-Markovian 
electron dynamics in a random potential ("memory ef- 
fects"). This classical contribution to the OPC (in con- 
trast to the classical contribution to the oscillatory ac 



conductivity) was shown 41 to be parametrically smaller 
than the quantum contributions we consider below. 

B. Quantum mechanisms of the oscillatory 
photoconductivity 

The Landau quantization leads to a periodic modula- 
tion of the DOS for high Landau levels, v(e) — v(s-\-uj c ), 
which, as far as the OPC is concerned, essentially modi- 
fies the classical picture above. In particular, in the high- 
temperature limit (T ^> (jj,u) c ), the distribution function 
acquires fast energy oscillations with the period u c on 
top of a smooth thermal smearing (see Appendix A): 

f(e,t,tp) = f T (e)+J r (s,t,ip)d £ f T {e), 
T(e,t,Lp) = f(£ + w c ,f,<?) , (3.3) 

where /r(e) describes the thermal distribution. What is 
important to us is that the amplitude of the part of T 
that oscillates with e oscillates also with the ratio oj/co c , 
which gives rise to the OPC. 

Neglecting the effect of microwaves on the DOS, we 
thus assume that the function gf+ Tt t (ip) does not depend 
on t and ip. Then, using Eq. (3.1), the dissipative dc 
current (2.18) can be expressed as 

-to,/**)^^), (3.4) 

where Fio(s) is the time- independent first angular har- 
monic of the distribution function (1.1). 
In the absence of external fields, 

Stt ph) {f(e)} = (3.5) 

and the 2DEG is in equilibrium with a thermal bath at 
temperature T, /(e) = /t(e)- Applying the electric field 
produces perturbative corrections to fx'- 

f = fT + ^2(C- 1 St im ) n f T , (3.6) 

n 

where C^ 1 = (dt + uj c d v — St; n ) _1 is the propagator in 
the l.h.s. of the QBE (2.11) and the collisions integral 
St; m , given by Eq. (2.13), should be expanded in powers 
of E dc and E w according to Eqs. (2.21)- (2.27): 

C + _ c,(a,m,k) 
Otim ~ im ' 

mk 

st (a, ro ,fc) a £ ( rn, k) ^ (g ?) 

The index a = _L, j (introduced for ease of visualization) 
shows which part of the kernel, Kj or K,^ , enters Stj m at 
a given order in the fields: a — _L or j for even and odd 
m + k, respectively. 
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The condition r q <C n r allows us to neglect K. v and 
higher-order terms in the expansion of the kernel (2.25). 
To the same accuracy, we omit all terms in Eq. (3.6) with 
ICj entering more than once, since these are small in the 



same parameter T q /r tr . The terms of zero order in K,j 
do not contribute to the current, as they are even in the 
electric fields and, hence, in (p. As a result, we represent 
-Fio, which determines the dc current (3.4), in the form 



F 



10 



E E 

vn M+A'=odd 



L E 

fi-\~X— even 



st 



c- 1 St 



U,M,K) 



1 E 

m+fc— even 



St 



(_L,m,fc) 



/t 



(3.8) 



where the summation indices run over nonnegative values 
and fi + x an d m + k are even, while M + K is odd. 

To obtain the linear dc response in the absence of radi- 
ation from Eq. (3.8), we put M — 1 and all other indices 
equal to zero, which yields 



F 



10 



(e-*(u B d v )-iSi& lfi) {f T }) ■ (3-9) 



Here St^ 1 ' ^ {fx} is given by Eq. (2.13) with IC repre- 
sented by ttj, Eq. (2.27), and the latter taken at first 
order in the dc field and at zero order in the microwave 
field, 



K 



(1,0) 



(TqT tr 



,1/2 



(3.10) 



Being substituted into Eq. (3.4), Fio from Eq. (3.9) gives 
the familiar linear-response dc conductivity, including 
the Shubnikov-de Haas oscillations and the nonoscilla- 
tory quantum correction [see Eqs. (4.14) and (5.7) be- 
low]. 

The effect of the microwave field on the dc current 
appears at order 0(Ed c E^). To this order, we pick up all 
terms with [iv + M + mn = 1 and \ v + K + kn = 2 in the 
expansion (3.8). The result is illustrated in Fig. 1, where 
diagrams (A)-(F) correspond to the following terms in 
Eq. (3.8): 

• "displacement" contribution (A): M = 1, K = 2, 
v = n = 0; 

• "inelastic" (B) and "quadrupole" (C) contribu- 
tions: M = n = 1, k = 2, K = v = m = 0. 
The two contributions are distinguished by the an- 
gular dependence of St[ r ^'°' 2 ' ) : (B) and (C) include 



(0,2) 



respec- 



the isotropic and quadrupole parts of K, 
tively; 

• "photovoltaic" contributions (D): K 
k = 1, M = v = 0; 



• diagrams (E): M = i/ = l,x = 2, K = n = n = 0; 

• diagrams (F): K = v = fx = x = ^-i M = n = 0. 
We now analyze these terms in more detail. 



(A) 



(0,0) (1,0) 




(0,0) (1,0) 



(C) 




(0,0) (1,0) (2,0) 




(2,1) 



(0,1) 




(0,0) (1,0) 



(0,0) (1,0) 



(E) 




+ 




(0,0) (1,0) 



(-1,0) (0,0) (1,0) 



(F) 



(1,1) 



(-1,1) 





(0,0) (1,0) 



(0,0) (1,0) 



FIG. 1: All possible contributions to the photoconductivity 
at minimal order O(EdcE^) are illustrated as graphs in (v,n) 
space of angular and temporal harmonics of the distribution 
function F vn , according to Eqs. (3.4) and (3.8). The arrows 
denote coupling of the harmonics F vn by the impurity colli- 
sion integral Sti m [Eq. (2.10)]. The filled circles denote the 
action of the propagator £ _1 in the resulting state. More pre- 
cisely, StJ^ 1 ' 2 ' (which is of third order in the electric fields) 
is denoted by the thick line in diagram (A), while the first- 
order collision integrals St;*i and StJ^ 0,1 ' in diagrams (B)- 
(F) are denoted by the thin lines. The double-line arrows 
in diagrams (B)-(F) correspond to the second-order integrals 



St 



(J-.M) 



and St 



(±,0,2) 



1. Displacement mechanism (A) 

The displacement contribution to the distribution 
function is obtained by expanding K,j to order 0(E t \ c E^ J ), 



(1,2) _ 



7 rrrt (9 v Xd c X u +Xd c X u d ip ) , (3.11) 
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and averaging X% over t 

-TTTT 1 



uit. 



X l = g ( £ + + £ - + 2£ + £ - cos2( ^) sin — ' ( 3 - 12 ) 

while neglecting all effects related to JC± by putting v = 
n = in Eq. (3.8). The result, 



(3.13) 



is represented in Fig. 1 by diagram (A). The thick arrow 
denotes the action of the term in JQj ' proportional to 
exp(i</?), which couples the zero angular harmonic /oo = 
It to /i . The filled circle stands for the propagator 
C^ 1 = (cjJcdtp)^ 1 in the resulting state (1,0). 

Equation (3.13) can be understood as follows. In 
the absence of the microwave field, the dc current oc- 
curs due to a difference in the rates of disorder-induced 
scattering with electron displacements along and against 
the applied dc field, so that the distribution function 
becomes anisotropic, Eq. (3.9). Absorption or emis- 
sion of microwave quanta during these scattering pro- 
cesses, Eq. (3.13), modifies the current by changing 
the final states and their occupancy. The displace- 
ment mechanism of the OPC was predicted long ago 
in Ref. 32 and addressed in the majority of theoretical 

works 32,33,36-39,44,45,47,48,55-57,61-65 ^ ^ Qpc ^ 

prehensive study is presented in Ref. 39. 



2. Inelastic mechanism (B) 

The inelastic contribution to the photoconductivity is 
governed by a change of the isotropic time independent 
part of the distribution function, i<oo , due to the absorp- 
tion and emission of microwave quanta. The change of 
F 00 is related to the t and tp independent part of the 
kernel, which, at order 0(E 2 ) 1 is given by 



X 2 , 



= "2^ 



c 2 \ -2 

£_ j sin — — 



(3.14) 



The QBE for F 00 to order 0(E%) reads 



St 



(e-ph) , 



{^oo} + St[r e) {^oo} = -(St 



,(-L,0,2) 



{/t} 



(3.15) 



Equation (3.15) with Eq. (3.14) substituted into the r.h.s. 
describes two effects: (i) electron heating, i.e., a mod- 
ification of the smooth part of the electron distribu- 
tion, and (ii) the appearance of an oscillatory correction 
to the isotropic part of the distribution on top of the 
smooth part [cf. Eq. (3.3)]. The heating is controlled by 
electron-phonon scattering [the first term on the l.h.s. 
of Eq. (3.15)] and leads to an effective electronic tem- 
perature T e larger than the bath temperature T. We 



relegate the discussion of T e as a function of the system 
parameters and the microwave power to Sec. VIII. Much 
faster processes of electron-electron inelastic scattering 
(the second term on the l.h.s.) cannot stabilize the in- 
crease of T e , but are capable of equilibrating electrons 
among themselves. The amplitude of the part of T os- 
cillating with e in Eq. (3.3) and the amplitude of the os- 
cillations in the photoconductivity are thus proportional 
to the rate of e-e collisions, 50 t^ 1 = T^(T e ), 



St 



(e-e) 



{^oo} 



00 



(3.16) 



Having obtained the nonequilibrium F o, we use ICj given 
by Eq. (3.10) to calculate the linear response with respect 
to the dc field, so that altogether at order 0{EfcE^) we 
have: 



r 10 — 



L St 



St 



(J-,0,2) 



{/t} 



(3.17) 

Equation (3.17) is represented by diagram (B) in 
Fig. 1, where the double- line loop denotes the action of 

( St^°' 2) ) on St and the thin arrow corresponds to the 

action of St-^ 1,0 ^ on the resulting Foo- 

One can see that, as compared to the displacement con- 
tribution to the OPC, the effect of the inelastic mecha- 
nism is accumulated during a much longer time r; n 3> r q . 
The period and the phase of the oscillations are the same 
for the two mechanisms (and correspond to those ob- 
served in the experiment); however, the amplitude of the 
inelastic contribution is r; n /T q >• 1 times larger. An- 
other important difference is that the T dependence of 
ri n makes it possible to explain the temperature-induced 
decay of the oscillations as observed in the experiment, 
in contrast to the oscillations produced by the displace- 
ment mechanism, which are T independent. The inelas- 
tic mechanism of the OPC was proposed in Ref. 34 and 
discussed in detail in Ref. 50. 

While mechanisms (A) and (B) are widely discussed 
in the literature, already at first order in the microwave 
power there exist additional contributions (C)-(F) to the 
OPC, which have not been studied before. 



3. Quadrupole mechanism ( C) 

In addition to the effect on Fqq, the absorption and 
emission of microwave quanta, described by iC^ 2 ^ - 

— — 1 



-r q 1 X 2 , leads to the appearance of a quadrupole cor- 
rection F20 to the distribution function [see Eq. (3.12)]: 

F 20 = (e- 2 ^( Wc ^)- 1 St™ y {M\ . (3.18) 



Acting by the operator A^ 1,0 ' [Eq. (3.10)] on the 
quadrupole correction yields 

F< C > - ( e-*^)- 1 St<* 1,0) {e 2 * J F 20 } ) , (3.19) 

as illustrated by diagram (C) in Fig. 1. 
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4- Photovoltaic mechanism (D) 

The photovoltaic contribution to the photoconductiv- 
ity is generated by the action of 



K 



(i,i) 



-2t 1 X UJ Xdc 



(3.20) 



on the equilibrium distribution, which leads to the ex- 
citation of the components F 01 and F 2 \ oscillating with 
time. Acting then by 

Jtf' 1] = -i(r q r tr )- 1/2 (9 v Af w + X„d v ) , (3.21) 

we calculate the ac response in the resulting state, which 
gives 



F< D > = (e-*^)" 1 



x St 



(j',o,i) 



(^ + 9 t )- 1 st i ( ^ M) {M 



(3.22) 



see diagrams (D) in Fig. 1 



5. "Inverse-order" contributions (E) and (F) 

Additionally, contributions similar to (B)-(D) but with 
the inverse order of the operators St^ and St^ in 
Eq. (3.8) are also possible. Diagrams (E) and (F) in 
Fig. 1, which describe these processes, involve harmon- 
ics of the distribution function different from those in 
diagrams (B)-(D). Namely, while contributions (B)-(D) 
include even angular harmonics only, except for the final 
Fio , all harmonics in diagrams (E) and (F) are odd in the 
angle of the momentum, except for the initial Fqo = Jt- 
Despite yielding nonzero contributions to F w , diagrams 
(E) and (F) vanish in the dc current in the limit of high 
T and strongly overlapping Landau levels [see the para- 
graph following Eq. (4.11)]. The explicit calculation in 
Sees. IV and V shows that, in the limit of strongly over- 
lapping Landau levels, mechanisms (A)-(D) give a com- 
plete set of contributions to the photoresponse to leading 
order in the microwave power. 



IV. QBE FOR OVERLAPPING LANDAU 
LEVELS 

From now on, we study the effect of microwave radia- 
tion on dc transport in the case when the Landau quan- 
tization modulates the DOS only weakly. In this limit, 
the modulation is represented by a single cosine term, 

i>(e) = v(e)jvQ = 1 — 2(5 cosets , t B — 2n/\u) c \ , (4.1) 

where 

5 = exp(-t B /2r q ) « 1 . (4.2) 



It is worth mentioning that the modulation of the DOS 
at order 0(5) is insensitive to the external electric fields 
[corrections 39,78 to the DOS induced by the electric fields 
appear at order 0{8 2 )}. 

The existence of the small parameter 8 greatly simpli- 
fies the solution of the kinetic equation in energy space, 
see Appendix B. The solution can be represented in the 
form (3.3) with 

F{e,t,ip) = <t> (<p, t) 

- 28 Re [</>i((p, t)exp(iet B )] ■ (4.3) 

For T 3> tg , the integration over e in Eq. (3.4) averages 
out terms in the current of first order in S. As a result, 
in the high-T limit the oscillations of the DOS manifest 
themselves in the current at order 0(S 2 ). We therefore 
expand 0o up to a term quadratic in 8: 



25- 



2^(2) 



(4.4) 



where 4>o is the Drude part independent of 8, and write 
the current to order (D(S 2 ) as 

Yd = 2ev F v ( e-* + 2S 2 Re 7^) ) . (4.5) 



It is convenient to introduce new functions 
according to 



and 



- 



4 2) +Re0i, 
Im0i , 



(4.6) 
(4.7) 



after which the current (4.5) is split into the classical 
(Drude) and quantum parts as follows: 



j d = 2a°E dc + 4S 2 ev F v (^j ( e* & ) 



(4.8) 



As shown in Appendix B, the functions 4>q , cfij, and <f>±_ 
obey 



i(d t + oj c d v )4>Q = d t _ JCj(t, i_) 



t_=o 



i(d t + uj c d^ 



(d t _ - -dt)Kj{t t t-) 



— fcj(t,tB)<f>± - <f>±£j{t,tB) , 

[d t +uj c d v - JC±(t,t B )](j)± +rr n 1 (q 



(-d t -d t _)ic ± (t,t- 



(4.9) 

'■B 

(4.10) 
(4.11) 



Here and below, any operator ICj t ± in the rightmost po- 
sition is understood as acting on unity. 

The first term on the r.h.s. of Eq. (4.10) gives the 
displacement contribution (A) to the dc current (4.8), 
while the second and third terms yield the inelastic, 
quadrupole, and photovoltaic contributions (B)-(D). It 
is important to notice that Eq. (4.10) for the function 
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<j>j includes the odd part of the kernel JCj only, while the 
function <fi± [Eq. (4.11)] is completely determined by the 
even part JC±. It follows that diagrams (E) and (F) in 
Fig. 1 produce no contribution to the dc current, because 
of the absence of coupling between <f>j and K.±. Although 

(2) 

the functions <p y and Re^i are coupled with each other 
[see Eqs. (B.6) and (B.8)] and do have contributions of 
type (E) and (F), each of these contributions vanishes 
in the combination [Eq. (4.6)] which determines the 
current according to Eq. (4.8). 

As follows from Eq. (4.8), for the calculation of we 
need only the t independent parts of and <j>j . Using 
the explicit form of ICj [Eq. (2.27)] and the fact that 
X(t- = 0) = 0, the r.h.s. of Eq. (4.9) averaged over t 
reads 



ft_ JCj(t,t-) 



-id,. 



d f X 



dc 



t_=0 



which gives 



To ev F f cos ip 

-*>dc 



sm ip 



and, correspondingly, the Drudc conductivity 



2 2 

2uj 2 T tI 



Similarly, 



where 



-ft B + 2</> ± X 



uJc^/t^ (1 + A 2 ) 



2 ' 



X = X\t_=t B — X dc + A a 
and for later use we introduce also 



A, 



X u , 



li_=t B 
Adc = <%dc| t _=t. 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



The derivative dt B in Eq. (4.15) should be understood 
according to the notation which we also use below in the 
paper: 

d tB A(X u ,X dc ) = d t _A(X u ,X dc )\ t _ =tB , (4.18) 

where A is an arbitrary function. 

Now that <pj, which determines the dc current (4.8), is 
obtained in the form of Eq. (4.15), the problem of finding 
the OPC to arbitrary order in E 2 or E dc is reduced to 
solving Eq. (4.11) for the function <f>±. In this paper, we 
restrict ourselves to the calculation of the linear response 
with respect to the dc field. In Sec. V, we obtain explicit 
expressions for the linear photoconductivity, correspond- 
ing to j oc E dc E 2 . With increasing microwave power, the 
system passes through several nonlinear regimes, which 
are studied in Sees. VI, VII, and VIII. 



V. LINEAR PHOTOCONDUCTIVITY 

According to the general classification presented in 
Sec. Ill, at order 0(E dc E 2 ) one can distinguish four dif- 
ferent mechanisms of the photoconductivity, illustrated 
by diagrams (A)-(D) in Fig. 1. We now calculate terms in 
the photoconductivity associated with the corresponding 
noncquilibrium corrections to the distribution function, 
given by Eqs. (3.8), (3.13), (3.17), (3.19), and (3.22). 

The displacement contribution (A), Eq. (3.13), is pro- 
duced by JCf a) cx E dc E 2 . All effects on even angular 
harmonics of the distribution function, governed by K.± , 
are neglected, i.e., <j)^ — 0. Equation (4.15) then re- 
duces to 



6ft 



AdcA2 . 



(5.1) 



By contrast, contributions (B)-(D) are related to the 
effect of the microwave and dc fields on even angular 
harmonics. At order 0(X 2 ), Eq. (4.11) reads 

(ft +u Jc d v )4>i_ +T,- 1 (<M - r-^ft, - ft/2)A 2 . (5.2) 

To linear order in the dc field, we replace X 2 in the 
r.h.s. of Eq. (5.2) by X 2 + 2X UJ X dc . The isotropic and 
quadrupole parts of X 2 [Eq. (3.12)] yield the inelastic (B) 
and quadrupole (C) contributions, respectively, while the 
photovoltaic contribution (D) is generated by the prod- 
uct A w A dc : 



^ = ^ dtB (X 2 ) 



b (c) 



d tB (— — 



U) c Tqd (p 



(xi - {XI)) 



if ) = T- 1 (d t + cj c d v )- 1 (2d tB - d t )x u x dc . 



(5.3) 
(5.4) 

(5.5) 



In Eqs. (5.3)-(5.5), 4>)_ ' and (f>\_ ' are averaged over t, 
whereas <j>^ is not. The current (4.8) is then produced 
by the dc [in the case of (B) and (C)] or ac [in the case 
of (D)] response in the resulting state: 



— ■ — 2 



ij_(A dc +X U ) . 



(5.6) 



Substituting Eqs. (5.3)-(5.5) into Eq. (5.6) and using the 

resulting f, together with <j>( A) [Eq. (5.1)] in Eq. (4.8) 

yields the current j = o- p hEd c - 

It is convenient to parametrize the photoconductivity 
tensor <r p h by four functions d s , d a , h s , h a : 



gph 

2cr D 



25 2 



1 + 25 2 -LU c T tr 

UJ c T tr 1 + 2S 2 

d s h a 

-h a d s 



28 2 



d a h s 
h s -d a 



(5.7) 
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The first matrix in Eq. (5.7) is the linear dc conductivity 
in the absence of microwaves. The quantum correction 
2<5 2 to the dissipative part comes from inserting Eq. (4.15) 
into Eq. (4.8) if one puts <p± = in Eq. (4.15) and also 
substitute X& c for X/(l + X 2 ) 2 . The resulting <f>j coin- 
cides with [Eq. (4.13)]. 

Two other matrices in Eq. (5.7) represent the 
microwave-induced corrections. We express the elements 
of these matrices in terms of the following dimensionless 
parameters: 



w = = ttlo/\lo c \ , 

Q = sin 2 u;(£ 2 +£ 2 ) , 

Qai = sin 2 w£ + £- , 

Qs = 2 sin 2 w 



2 LU +UJ C 2 CO - U> c 



LO + 2lO c 



lo — 2lo c 



Qas 



2 sin z w — 

LO 

2 1 



LO + 2lo c 



E 2 



LO 



2co r 



(5.8) 
(5.9) 
(5.10) 

(5.11) 



(5.12) 



where £± are defined in Eq. (2.24). In addition to the 
effective microwave power Q we introduced three sim- 
ilar parameters, Eqs. (5.10)-(5.12), with the subscripts 
"AI","S", and "AS" standing for "anisotropic", "sym- 
metric", and "antisymmetric", respectively (in accord 
with their appearance in the corresponding components 
of <7 p h, see below). 

The diagonal isotropic part d s of <r p h, Eq. (5.7), is a 
sum of the displacement, inelastic, and photovoltaic con- 
tributions: 



d( A) = 3Q + 6Q w cot w , 
d (B) = 2 ^Q WC oti 



;w 



d{^ = 



wQ s 



(5.13) 
(5.14) 
(5.15) 

(5.16) 



The diagonal anisotropic part d a is governed by the dis- 
placement and photovoltaic mechanisms: 

da = 4 A) + 4 0) , 

d[ A) = 6Q A i(l + 2wcotw) , 



d^ 



4w Qai 



UJT a 



(5.17) 
(5.18) 

(5.19) 



The nondiagonal symmetric term h s is due to the 
quadrupole mechanism: 



h s = = ^ AI w cot w 
w c r q 



(5.20) 



while the antisymmetric Hall part is generated by the 
photovoltaic mechanism: 

ha = /ii D) = — (1 + wcotw;) . (5.21) 



Let us now discuss the results (5.13)— (5.21), obtained 
for the linear regime in the microwave power: 

• the inelastic mechanism yields the dominant effect 
because of r in /T q > 1: (B) > (A), (C), (D); 

• (C), (D) ^> (A) for the case of cj c r q ~ ujr q <C 1 
(strongly overlapping Landau levels). Note, how- 
ever, that since 5 [Eq. (4.2)] depends exponentially 
on lo c , a pronounced effect (from the practical point 
of view) is only possible at not too small magnetic 
fields. The effects of (A), (C), and (D) are then 
comparable in magnitude, differing only in In S; 

• although (B)»(C),(D), the oscillatory correc- 
tions to the Hall conductivity are entirely governed 
by the quadrupole (C) and photovoltaic (D) mech- 
anisms; 

• the quadrupole contribution h s violates Onsager 
symmetry: <J ik {B, s_, s+) ^ au{-B,s + ,S-) [the 
effect of time inversion on the microwave field po- 
larization can be expressed as s + s_, which 
is equivalent to to — > -co, see Eq. (2.24)]. Al- 
though Onsager symmetry need not be satisfied in 
nonequilibrium systems, there exist not many ex- 
amples of its violation. One of them, for a dif- 
ferent system but under similar conditions, was 
found in Ref. 77. Symmetry with respect to the 
parity transformation, x —* —x, must be satis- 
fied: a ik (B,S-,s + ) = r ik a ik {-B,s + ,S-), where 
^xx — ^*yy — L ^*%y — ^y% ^* ^ can be eas- 
ily verified in our case since all parameters (5.8)- 
(5.12) are unaffected by the simultaneous change 
B — > —B and s + s_; 

• the quadrupole term h s and the anisotropic term in 
the diagonal photoresponse d a vanish in the case of 
circular polarization of the microwaves (Qai = 0) . 



VI. FEEDBACK EFFECT: SATURATION OF 
THE INELASTIC CONTRIBUTION (SIC) 

The results for <r p h oc E 2 obtained in Sec. V, are valid 
at not too high microwave power. The QBE in the form 
(4.11) contains two sources of nonlinearity in E 2 . At 
Q > 1, high-order terms in the expansion of ftj and JC± in 
powers of E u become important, leading to multiphoton 
corrections to both </>j_ and tfij, see Eqs. (4.11) and (4.15). 
The other source of nonlinearity, the "feedback" term 
JC±(f)± in Eq. (4.11), strongly modifies the photoresponse 
at much smaller microwave power, Q m ~ 1, where 



Qin — „ Q • 

2t q 



(6.1) 



Since T; n S> r q , the multiphoton corrections can still be 
neglected at Q m ~ 1. As we show below in this section, 
at Qi n » 1 the single-photon feedback effect leads to a 
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saturation of the inelastic contribution and to a strong 
modification of all other contributions to the photocon- 
ductivity. 

A. Feedback effect on the isotropic part of the 
distribution 



phonon bath: the dynamic equilibrium between these two 
processes determines the effective electronic temperature 
T e , see Sec. IX. However, the saturation (6.5) of the 
amplitude of the oscillations in the isotropic part of the 
distribution function is governed by the same feedback 
effect as the transparency in the two-level system. 



The feedback effect on the inelastic mechanism was 
studied in detail in Ref. 50. Here, we reproduce the 
results 50 within the more general framework of Sees. Ill 
and IV. Using Eq. (2.26), we rewrite Eq. (4.11) in the 
limit Q <C 1 as 



0J 



.) 2d, 



Lb 



di 



2r n 



X 1 



(6.2) 



Averaging Eq. (6.2) over t and 4>, we obtain the equation 

(B) 

for the isotropic time-independent part 4>\ = (<j>±): 



b (B) 



= d tl 



(6.3) 



In the linear regime with respect to the microwave power, 
i.e., for Q <C T q /n n , we return to (f>\_ oc Q given by 
Eq. (5.3). At Q <~ T q /r in <C 1 two terms on the l.h.s. 
of Eq. (6.3) become of the same order of magnitude, 
whereas the last term on the r.h.s. is still small and 
can be neglected (corrections of higher order in Q, gen- 
erated by this term, are considered in Sec. VIC). We 
thus obtain for Q <C 1 and arbitrary order in Q ln : 



i + Oi 



-to cot w 



(6.4) 



It follows that in the limit Q m > 1 we reach a saturation 
of the inelastic contribution (SIC): 



(6.5) 



The dc response in the stationary nonequilibrium state 
(6.4) yields 



4Qi, 



i + Oi 



■ w cot w 



(6.6) 



reproducing Eq. (15) of Ref. 50. 

The saturation of the inelastic contribution at large 
Qi n has much in common with the well-known effect of 
"self-induced transparency" for a two-level system in a 
strong resonant electromagnetic field, which occurs due 
to equalization of the population of the levels. Clearly, in 
the system we consider here the complete transparency 
cannot be obtained since the electron spectrum in our 
problem is continuous. Even in the SIC regime, absorp- 
tion does not tend to zero. Electrons still absorb mi- 
crowave quanta and transmit the absorbed power to the 



B. Quadrupole, photovoltaic, and displacement 
contributions in the SIC regime 

In the nonlinear regime, the strong oscillations of 
Fooi?) produced by the inelastic mechanism essentially 
modify contributions (A), (C), and (D). Since (f>)_ and 
<t>^ are still small at Q m <~ 1, we replace (j>± in the 
term X 2 cf)j_/Tq on the l.h.s. of Eq. (6.2) by (f>^ . Solving 
Eq. (6.2) for the quadrupole and photovoltaic contribu- 
tions then gives: 



b (DB) 



LU c Tq \ / 

2{d t + Lo c d v )- 1 

d tB - d t /2 



XujXdc 



(6.8) 



Here and below we emphasize the influence of the in- 
elastic mechanism by adding the second superscript B. 
Substituting Eqs._(6.7) and (6.8) into Eq. (4.15) and us- 
ing the resulting <j>j in Eq. (4.8) yields the photovoltaic 
and quadrupole contributions to the Hall part of <j p h 
[Eq. (5.7)]: 



h (CB) _ Qai wcotw 

W C T q Qin + 1 ' 
^{DB) _ QaS A Qin 



1 



Qin + 1 



-w cot w 



(6.9) 



(6.10) 



The diagonal photovoltaic contributions di D ^ and d!^\ 
which originate from the derivative d t in Eq. (6.8), are 
not affected by the oscillations in F o and remain the 
same as in the linear regime, Eqs. (5.16) and (5.19). 

The oscillations of F o also modify the displacement 
contribution to the current. Similarly to Eqs. (6.7) and 

(6.8) we substitute for cf>j_ in Eq. (4.15), which gives 



Using Eq. (6.11) in Eq. (4.8) we have 



&x Ac xi . 



= 3Q ( 1 



1 



Qin + f 



2w cot w 



(6.11) 



(6.12) 



(6.13) 
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We thus see that in the SIC regime, Q m 3> 1, the 
quadrupole term lif 3 ^ saturates as a function of the 
microwave power, whereas the nondiagonal photovoltaic 
term h^ 3 ^ and the both displacement terms di AB ^ and 
da 45 ' 1 continue to grow linearly with increasing power but 
are strongly modified as compared to the linear regime. 



C B) 

up the contributions (6.14) and (6.15) to <pj_ we obtain 
di B ^ — 2tB(f>^ with the linear-in-Q terms included: 



4w Qjr 



cotw 1 



Q 2 



3 g Q in + 2 

4 Qin + 1 / 



q 

(6.18) 



C. Two— photon and double— frequency corrections 
to the isotropic part of the distribution 

For Q in » 1, the main (inelastic) contribution to the 
oscillatory part of the photoconductivity is independent 
of the microwave power as long as Q is sufficiently small. 
Deviations from the "plateau" come from contributions 
(A) and (D) (which, as shown in Sec. VI B, grow linearly 
with Q in the SIC regime), and also from corrections of 

order O(Q) to the Q independent term cjy^ [Eq. (6.5)], 
which we obtain below. 

Equation (4.11) produces two types of corrections to 

4>\ ' of linear order in Q. One of them is related to two- 
photon processes in K,±, which come from the expansion 
of IC± to fourth order in X u , yielding 



W + {XI Xt) 



= U) cot w ■ 



1 _^ Q Qm + 2^ (614) 



Qin + 1 V 4"Q in + l 

The other is associated with single-photon excitation of 
the second temporal harmonic of the distribution func- 
tion. This correction originates from the term dtX 2 in 
Eq. (6.2) (which has played no role so far) and to leading 
order in X 2 , reads: 



(6.15) 



where 



X-l-Xl = -sm 2 w[£lcos(2u;t + 2ip) 

+ £ 2 cos(2wt - 2tp) + 2£ + £_ cos 2ut ] . 

(6.16) 

The constant of integration C in Eq. (6.15) should be 
found from the periodic boundary condition for the func- 
tion <j)^\(p,t). Averaging Eq. (6.2) over ip and t, we 

get the boundary condition in the form T ln (X 2 (f>^^) — 

~ T q( ( f>± UJ) }, which gives C = T~ 1 Q 2 uQin/ (l+Qin), where 



Q 



sin 4 w 



2u> 



8Q \ uj + u c T uj — uj c 



£ 4 _ 



Equation (6.18) tells us that for strongly overlapping 
Landau levels the double-frequency correction is much 
stronger than the one coming from the two-photon pro- 
cesses, except for the close vicinity of the cyclotron- 
resonance harmonics where the two-photon correction ex- 
hibits singular behavior. 



D. Discussion of the results 

Equations (6.9), (6.10), (6.12), (6.13), and (6.18), to- 
gether with Eqs. (5.16) and (5.19), describe the photo- 
conductivity <7ph [Eq. (5.7)] to arbitrary order in Qin = 
T- m Q/2T q and to first order in Q. In the limit Q in <C 1, 
Eqs. (5.13)— (5.21) are reproduced. At Qin ^ 1, the evo- 
lution of the photoresponse with increasing power can be 
summarized as following: 

( B) 

• Growth of the leading (inelastic) term d s satu- 
rates due to the feedback effect, Eq. (6.6), and so 
does the quadrupole term hi C3 \ Eq. (6.9). 

• Remarkably the parts proportional to wcotw of 
the displacement contributions di AB ^ and d^ B \ as 
well as that of the photovoltaic contribution h\F B \ 
change sign at Qi n = 1, and the linear growth is 
recovered with opposite sign at higher power. 

• While at Qi n -C 1 the inelastic contribution is the 
only one that is T dependent, in the crossover re- 
gion, Qi n <~ 1, all contributions, except for the di- 
agonal photovoltaic terms d{ D ^ and di°\ change 
strongly as a function of T. 

• However, apart from the small quadrupole cor- 
rection h\ , the photoconductivity in the limit 
Qin S> 1 becomes independent of T; n and, hence, of 
temperature: 



(6.17) 



d^ 


= w cot to (4 - 3Q) + 4wQ 2a) /wT q , 


(6.19) 


JAB) 


= 3Q (1 — 2wcotw) , 


(6.20) 


d iAB) 


= 6Qai(1 — 2wcotw) , 


(6.21) 


d^ 


= wQ S /u}T q , 


(6.22) 




= -4KlQ A l/wTq , 


(6.23) 


h iCB) 


= (QAl/Qin)™ COt w/uJ c T q , 


(6.24) 


h V>B) 


= Qas(1 - w cot w)/uj c T q . 


(6.25) 



Only the constant C contributes to the double- frequency 
correction to the dc current at first order in Q. Summing 



Equations (6.19)-(6.25) predict that there are two 
sources of nonlinearities which will essentially modify 
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the plateau behavior in the SIC regime with increas- 
ing microwave power: the multiphoton processes and 
the "spreading" of the feedback effect over higher tem- 
poral and angular harmonics. The multiphoton effects 
become strong at Q > 1, whereas the feedback- related 
excitation of distant harmonics becomes strong already 
at much lower power, namely at Q > uj c r q . It is impor- 
tant that the strong deviations from the plateau behav- 
ior take place when Q is already much larger than r q /T; n 
and, therefore, inelastic scattering can be completely ne- 
glected. In Sec. VII we discuss the regime of ultrahigh 
power in more detail. 



VII. ULTRAHIGH POWER: FEEDBACK VS 
MULTIPHOTON EFFECTS 

A. Classification of contributions to <r p h at 
ultrahigh power 

In the SIC regime considered in Sec. VI, strong oscil- 
lations in isotropic part of the distribution function es- 
sentially modified all contributions to the dc current (A)- 
(D). At ultrahigh power, contributions involving high an- 
gular and temporal harmonics of the distribution func- 
tion become important, so that our previous (A)-(D) 
classification, developed in Sec. Ill for the linear <j p h, 
should be generalized. Despite the excitation of the high 
harmonics, at ultrahigh power it is still possible to dis- 
tinguish the generalized displacement (A), photovoltaic 
(D), and non-photovoltaic contributions. In the linear 
regime with respect to the dc field E dc , the latter enters 
the part of the kernel JC± in the photovoltaic term only. 
In all other terms, E dc belongs to JCj. Further, in the dis- 
placement term all effects of the microwave and dc fields 
on even angular harmonics are absent. We now proceed 
to calculate various contributions to cr p h in the limit of 
ultrahigh power, where the inelastic scattering is of no 
importance. 

The displacement contribution is readily calculated to 
arbitrary order in microwave power by putting tf)j_ = in 
Eq. (4.15) for <f>j and extracting the part linear in E dc : 



evF 



Re (E+, 



x d tB t B 



1-3*2 
(1 + A2)3 



(7.1) 



where E^ c — (Ed c ):r + i(Edc)y and we used the explicit 
form of Adc given by Eq. (2.22). As for the effect of the 
ac and dc fields on even angular harmonics, Eq. (4.15) 
yields: 



2<p { f>X dc 



1 - 3X 2 



X u , 



p[D) _ 20 ± 



(7.2) 
(7.3) 



Here <f/^ is the solution of Eq. (4.11) that is independent 

of E dc (at small power, <(r®' is the sum of the inelastic and 

quadrupole terms) and (j/f* is the linear-in-i^dc (photo- 
voltaic) solution: 



d t + uJ c d v + 



1 X 2 



Tvi + xir 1 - 



a(°) 



d tB -d t /2 Xl 
r q 1 + X I 

1 xl 



_ 8 tB - dt/2 - «^ 0) 2X dc X L 



(7.4) 



(7.5) 



In Eq. (7.4), we omitted the term which describes in- 
elastic scattering, since in the following we consider the 
limit Q ^> T q /r in . In this limit, the inelastic scattering 
is irrelevant, as demonstrated in Sec. VI D. For simplic- 
ity, we focus here on the case of circularly polarized mi- 
crowaves. Linear polarization is briefly discussed in the 
end of Sec. VII C. 



B. Circular polarization: General solution 

For a circularly polarized microwave field, solution of 
Eqs. (7.4) and (7.5) greatly simplifies. For definiteness, 
let us consider passive circular polarization by putting 
s+ = 1, s_ = in Eq. (2.21), so that 



*w = \j~Q cos x , x = tp + uit , 



Q = — 

Ttr 



eE^VF 



LU(UJ C + LU) 



sin w . 



(7.6) 
(7.7) 



It is convenient to introduce new functions go{x) and 
gi(x) by casting <j)^(x) and <p^\x, if) in the form 



2i B 1 S , o(aO , 



evF 



x 2 Re 




-^dc 9i(x) exp(-iip) 



(7.8) 



(7.9) 



The partial differential equations (7.4) and (7.5) trans- 
form then into ordinary differential equations for go i i(x): 

g' N (x) + [ap(x) - i(3S N1 ] g N (x) = all^a;) , (7.10) 

where N = or 1, the constants a and (3 are given by 

a = 1/{uj + oj c )T q , (3 = uj c r q a , (7-11) 

and the functions p(x) and II jv (x) are written in terms 
of X u as 



p(x) 



V 



XI ' 



(7.12) 
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X 2 



n (x) = w(cot w + tanx) — ' (7-13) 

n 1 (x) = [i-2 9o (x)} {l *x l)2 

1 - 3X 2 



+ w{ cot w + tan x )X U _ 



(7.14) 

Integration of Eq. (7.10) with a periodic boundary con- 
dition #at(0) = <7jv(27r) yields 



9n{x)^Qn{x) + exp[-S N (x)] 



Gn(2tt) 



l-exp[-Sjv(27r) 



(7.15) 



X 

Q N (x) = a J dyU N (y)exp[S N (y) - S N (x)] , 
o 

(7.16) 

X 

S N (x) = -if3x5 N1 +a Jdyp(y) . (7.17) 
o 

Note that for the case of circular polarization the 
anisotropic contributions to the photoconductivity ten- 
sor <T ph [Eq. (5.7)] vanish, d a — h s = 0, so that <r p h 
reduces to 

. o d (l + 25 2 (l-d s ) -uj cTtI -2S 2 h a \ 

ph ~ xx \ oj cTtI + 25 2 h a 1 + 25\l -d,)J • 

(7.18) 

By using Eqs. (4.8), (7.1)-(7.3), (7.8), and (7.9), the 
terms d s and h a in Eq. (7.18) are expressed as 

d a = dW , (7.19) 

where 

1 - 3A 2 



d { s A) = l-(l + 2wcotwQd Q )( — 



df) =4Re ( gi ( 



(1 + *2)V, 



and 



/i a = -41m { 3i (x) 



(i + x 2 ) 2 



(7.20) 
(7.21) 

(7.22) 
(7.23) 



^ + X 2 ) 2 /x 
The brackets (. . .) x denote averaging over x. 

C. Photoresponse at ultrahigh power 

Being represented in the form of Eq. (7.20), the dis- 
placement contribution is easily calculated, 

d (A) = ! + r Q ' 2 -}^ + 2w cot w Q *~®{t , (7.24) 



reproducing Eq. (6.8) of Rcf. 39. Unlike di A \ the terms 
e?i°\ di D \ and h a , which are expressed through the func- 
tions go(x) and gi(x), cannot be written in a closed ana- 
lytical form for arbitrary a, (3, and Q. When discussing 
<Tph at ultrahigh power in various limiting cases, it is use- 
ful to represent d s and h a as 

d s = wcotwFid(Q,a,/3) 

+ w F 2d (Q, a, (3) + F 3d (Q, a, 0) , (7.25) 

h a = w cotwFih(Q, a, (3) 

+ w F 2h (Q, a, 13) + F 3h (Q, a, f3) , (7.26) 

where the different terms are in one-to-one correspon- 
dence with the terms in the functions IIv(x) on the r.h.s. 
of Eq. (7.10) proportional to w cot w, w, and independent 
of w (except for the dependence on lo in X u ). 

The results of a numerical calculation of d s and h a 
[Eqs. (7.20)-(7.23)] are illustrated in Figs. 2-5, where F nd 
and F n h (with n = 1, 2, 3) are drawn as a function of Q 
for (3 — 0.3 and several values of a = 1, 10, 30. Since Ti n 
is assumed to be infinite throughout Sec. VII, the scale 
of Q ~ T q /-ri n , below which the system exhibits the linear 
response with respect to the microwave power, vanishes. 
As a result, the function Fid is nonzero (Fid = 4) at 
Q = in Fig. 4. With the exception of F 3 d, ah other 
functions F n d and F n h show nonmonotonic behavior with 
increasing Q and decay to zero in the limit of large Q. 
The function F$d — > 1 for Q — > oo, thus eliminating, as 
can be seen from Eq. (7.18), the quantum correction 2S 2 
to the dissipative part of the classical Drude conductivity. 
It follows that all quantum effects in cr p h are destroyed 
and the classical Drude conductivity is restored in the 
limit of high power. 

For a>l (which is the case for strongly overlapping 
Landau levels, unless lo /lo c is too large), the displacement 
contribution (7.24) is clearly insufficient to correctly de- 
scribe the behavior of d s as a function of Q. This can 
be seen, e.g., in Fig. 2, where the displacement contribu- 
tion to Fsd is shown by the dashed line. A much faster 
behavior of F 3 d at small Q is due to the photovoltaic con- 
tribution. Indeed, the expansion of F 3d in powers of Q 
reads: 



F 3d ~ 3Q + 



(aQ) 2 1 + 3/3 2 
2 (1-/3 2 ) 2 



1 



a 



(7.27) 



where the linear term, due to the displacement mecha- 
nism, crosses over at Q <~ a~ 2 into the faster quadratic 
term related to the photovoltaic mechanism. At Q <~ 
ot , F^d reaches a maximum value of order unity. In the 
limit of large Q, F$d approaches unity from above. 

The displacement contribution to the function F 2d , 
shown in Fig. 3(a), is exactly zero, so that F 2d is entirely 
due to the contribution of d^ and d{ D \ At Q <C a -1 , 
F 2 d behaves linearly with Q: 



(1 + Q) 5 / 2 



(1 + Q) 7 / 2 



F 2 d ^ 



aQ 2aQ 



l-(3 2 



(7.28) 
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FIG. 2: Evolution of the contribution F 3d (Q), Eq. (7.25), 
to the photoconductivity (7.18) with increasing effective mi- 
crowave power Q, Eq. (7.7), for a = 1, 10, 30, and f3 = 0.3 
[Eq. (7.11)]. Dashed curve: the displacement contribution to 
F 3d (Q), given by the first two terms in Eq. (7.24). 



The first term in Eq. (7.28) comes from d^ and is asso- 
ciated with excitation of the double-frequency harmonics 
F 22 and F_ 2 ,-2 [Eqs. (6.15) and (6.19)]. The second, 
photovoltaic term, which is due to excitation of F 2 i, Fqi, 
F—2,-i, and Fo t -x, has a larger slope with opposite sign 
[Eq. (6.22)]. Their behavior at arbitrary Q is illustrated 
in Fig. 3(b): one can see that each term decreases slowly 
with increasing Q for Q 3> a -1 , but they strongly com- 
pensate each other, leading to a much faster decay of 
F 2 d- 

The function Fid is shown in Fig 4(a). It is finite at 
Q = and therefore is much larger than both F 2 d and 
at Q -C a -1 . The expansion of Fid in powers of Q 
reads: 



Fi d ~4-9Q 



{aQf 1 + 3/3 2 
2 (1-/3 2 ) 2 



(7.29) 



The nonzero value of Fid at Q = is due to the saturated 
inelastic contribution d s [Eq. (6.19)], whereas the linear 
term comes from both the two-photon correction to the 
inelastic contribution [Eq. (6.19)] and from the displace- 
ment contribution [Eq. (6.20)]. The two contributions 
related to and di A ^ are presented for arbitrary Q in 
Fig. 4(b). Similarly to F 3 d, the photovoltaic term grows 
quadratically at Q <C a -1 . As seen from Fig. 4(b), for 
a S> 1 it dominates at large Q, so that Fid changes sign, 
reaches negative values of order unity, and approaches 
zero in the limit Q — -> oo from below. 

The functions F n h contributing to the Hall part of <T p h 
are shown in Fig. 5. All three are related to the pho- 
tovoltaic mechanism and come from the imaginary part 
of gx(x) in Eq. (7.23). At Q < or 1 they grow with 
increasing microwave power, 



F2/1 



2/3 



1-/3 2 



(1-/3 2 ) 



aQ , 



(7.30) 
(7.31) 




FIG. 3: (a) Contribution F 2d (Q), Eq. (7.25), to the photo- 
conductivity (7.18) for a = 1, 10, 30, and (3 = 0.3; (b) "Non- 
photovoltaic" (positive) and photovoltaic (negative) contribu- 
tions to F2d{Q) for a = 1, 10. 



3/1 



1-/3- 



: aQ 



(7.32) 



and decay to zero in the limit of large Q. The maxi- 
mum values of |F ra h|, which are reached at Q ~ a" 1 , are 
of order unity (for /3 ~ 1), similarly to the behavior of 
F„d- It follows that at Q ~ a -1 the Hall part h a of the 
photoresponse becomes as strong as the diagonal part d s . 

The behavior of d s and h a in the limit of large Q can be 
found by observing that the main contribution to the in- 
tegrals (7.20)-(7.23) over x comes then from a close vicin- 
ity (of width Q^ 1 / 2 ) of points x = tt/2 and x = 3tt/2. 
Around these points the multiphoton processes, which 
sum up to produce the factors (l + X 2 ) -1 in Eqs. (7.20)- 
(7.23), are strongly suppressed. For Q 3> a 2 , the expan- 
sion of the factors exp(±Sjv) around unity in Eq. (7.16) 
yields the asymptotic behavior of the functions F„d and 

Fnh' 



Fid 

F3d 



la 1 



I6Q 3 / 2 ' 

™2 



8Q 3 / 2 
3/?a 
"4Q37 2 ' 

f3a 



2Q 3 / 2 



(7.33) 
(7.34) 
(7.35) 
(7.36) 



Note that the displacement mechanism gives terms of 
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(a) 
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2 



: = 10 



a = 30 



FIG. 4: (a) Contribution F ld (Q), Eq. (7.25), to the photo- 
conductivity (7.18) for a = 1, 10, 30, and /3 = 0.3; (b) Con- 
tributions to F ld (Q) [Eq. (7.19)] for a = 1, 10: Dashed line: 
Displacement contribution [third term in Eq. (7.24)]; solid line 
in the upper half-plane: nonphotovoltaic contribution [second 
term in Eq. (7.19)]; solid line in the lower half-plane: photo- 
voltaic contribution [third term in Eq. (7.19)]. 



order Q 3 / 2 in Fid and F%d, which are much smaller than 
the terms of order a 

g-3/2 in Eqg (j 33) and (j 34) for 

a > 1. 

In the intermediate interval a" 1 C Q C a 2 , where all 
the functions F n d and F n h, n — 1,2,3, fall off with in- 
creasing Q as a power law, the function gN(x) can be ap- 
proximated as gAr(x) ~ Gn{x) — TIn(x)/p(x) everywhere 
except for a vicinity of the singular points x = ir/2 and 
x = 37r/2. The singular behavior around these points 
requires special care for the function F2d, contributions 
to which coming from and only weakly depend 
on Q and strongly compensate each other, as illustrated 
in Fig. 3(b). Consider the part F^ related to d^: 



2d ir J w (1 + Qcos 2 x) 



1 — 3Q cos 2 x 



2 ~^3 



(7.37) 



where 



C(x) = aQ / dy 



sin y cos y 



(1 + Qcos 2 ?/) 2 
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FIG. 5: Hall contributions (7.26) to the photoconductivity 
(7.18) for a = I, 10, 30, and /3 = 0.3: (a) F lh (Q); (b) F 2h {Q); 
(c) F 3h (Q). 



x exp I — / dz 



aQ cos 2 z 
1 + Q cos 2 j 



(7.38) 



In Eq. (7.37) we neglected the exponentially small (for 
Q ^> a -1 ) boundary term in the solution (7.15). For 
Q ^> a^ 1 the main contribution to F^ given by 
Eq. (7.37) comes from the singularities at x = ir/2, 3ir/2. 
For Q«a ! one can neglect Q in all the denominators 

in Eqs. (7.37) and (7.38), after which F^ is represented 
as 



? (0) 

2d 



2CX Q f A f A ( \ 

— / dx\ I dx2 x\ sgn(aq — x-i) 



x exp — a 



- vs ■ (7M) 

The result does not depend on Q, which explains the 



17 



"plateau" in the behavior ol > (Q) shown in Fig. 3(b). 
It is important, however, that a similar calculation of the 
contribution to F 2 d coming from the photovoltaic term 
di D ^ exactly cancels the constant term in Eq. (7.39), so 
that the difference between the two falls off rapidly with 
increasing Q for Q » or 1 , as observed in Fig. 3(a). 

To summarize, in the limit of strongly overlapping Lan- 
dau levels and low temperatures, 1 <C a <C Ti n /r q , there 
appear five regions of Q with essentially different behav- 
ior of the photoconductivity as a function of Q and T: 

• Q <C T q /Tj n : the photoresponse is linear in Q and 
is strongly dependent on T; 

• Tq/^in <C Q -C 1/a: the strong feedback effect leads 
to the saturation of the photoresponse as a function 
of Q. In this interval of Q, as well as at higher Q, 
the photoconductivity ceases to depend on T; 

i a -1 « Q « 1: the feedback effects strongly sup- 
press the photoresponse; 

• 1 « Q « a 2 : the multiphoton effects become im- 
portant and modify the feedback effects; 

• a 2 <C Q: the photoresponse is dominated by the 
multiphoton effects, the feedback effects can be 
treated perturbatively. In the limit Q — > oo the 
classical Drude conductivity is restored. 

Above, in the bulk of Sec. VII C, we focused on the case 
of circular polarization in the limit of high power. Two 
main differences of linear polarization as compared to cir- 
cular polarization are as follows. First, the photoresponse 
becomes anisotropic, i.e., depends on the mutual orien- 
tation of E w and Ed c - Second, the photoresponse may 
show much stronger resonant features at noj = 2muj Cl 
where m and n are integer numbers. The latter is re- 
lated to the different structure of perturbation in (u, n) 
space [Eq. (1.1)] induced by the radiation. For the case 
of passive circular polarization considered above, only 
the harmonics f2m,2m along the diagonal are excited in 
the absence of the dc field [they are all contained in the 
function g in Eq. (7.8)]. The linear response to the dc 
field, included in Stj_, couples these harmonics with their 
neighbors i*2ro,2ro±i- Under the action of Stj all these 
harmonics contribute to F±i,o and thus to the dc photo- 
conductivity. In the case of linear polarization, however, 
all harmonics i^m.n, with m and n integer, are excited 
by radiation, which leads to the "spreading" of the per- 
turbation all over the (is, n) plane. The resonances at 
nuj = 2muj c for the case of linear polarization warrant 
additional study. 



VIII. MICROWAVE-INDUCED 
MAGNETORESISTIVITY OSCILLATIONS 

In Sees. V-VII we have calculated the photoconductiv- 
ity tensor <r p h and discussed its evolution with increasing 



"effective" microwave power Q defined in Eq. (5.9) as 



Q = P sin 2 — lo 2 



where 



(u> c + u) 2 (lu c -lu) 2 



2t„ / eE^vp 



7"t r \ OJ 



, (8.1) 



(8.2) 



However, what is usually measured in the experiment is 
the photoresistivity p p h as a function of uj c at fixed u) 
and P. The dependence of /5 p h on lo c /lo for a given P is 
discussed below, assuming, as in Sees. VII B and VII C, 
passive circular polarization of the microwave field. 

For circular polarization, the photoresistivity tensor is 
easily obtained in the limit of classically strong magnetic 
field, w c r tr ^> 1, by inverting Eq. (7.18): 



2 CT D> c T tr )2 

l + 25 2 (l-d s ) Lu c T tI - 26 2 h a 
-u; c r tr + 26 2 h a l + 25 2 (l-d s ] 



(8.3) 



Here we neglect a small admixture of — 46 2 h a /u> c Tt r in 
the diagonal part of p p h, as well as a similar admixture 
of ±4:S 2 d s /uj c Tt r in the Hall part. 

It was shown in Sec. VII that, while in the linear 
regime with respect to the microwave power the inelastic 
contribution to the photoresponse dominates, already at 



Q 



1 the inelastic and photovoltaic contribu- 



tions become comparable in magnitude. The displace- 
ment contribution becomes also relevant at this power 
if a ~ 1. The interplay of the three contributions in 
the traces of d s and h a as a function of the ratio u> c /lu is 
illustrated in Figs. 6-9 for the case of passive circular po- 
larization. The dependence of d 8 and h a on oj c is shown 
over a single period of oscillations around the second har- 
monic of the cyclotron resonance for 2/5 < u) c /oj < 2/3. 
Each of Figs. 6-9 corresponds to one of four values of the 
microwave power P — 0.01, 0.1, 1, and 10. 

Within each of Figs. 6-9 two cases are illustrated, cor- 
responding to two different values of the parameter a 
[Eq. (7.11)] taken at oj c /uj = 1/2. In every figure, (a) 
and (b) correspond to a.\ii = (3oj c T q ) _1 = 1, while (c) 
and (d) to a.\ji = 10. In all the figures we take the same 
ratio Ti n /r q = 50, thus neglecting the heating effect, dis- 
cussed separately in Sec. IX B. 

At P = 0.01 (Fig. 6), not only Q but also Q in = 
Ti n Q/2r q is small compared to unity in the whole interval 
2/5 < llJc/lj < 2/3. Therefore, the photoresponse is well 
described by the linear-in-P asymptotes (5.13)-(5.21). 
More specifically, the diagonal part at a ^ 1 is dominated 
by the inelastic contribution d s [Eq. (5.15)], while at 
a»l the photovoltaic contribution di D ^ [Eq. (5.16)] be- 
comes pronounced near the ends of the interval (where 
inelastic contribution vanishes, cotu> = 0). The Hall re- 
sponse is governed by the photovoltaic contribution ha 
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FIG. 6: Diagonal, d 3 , and Hall, h a , parts of the magnetore- 
sistivity tensor (8.3) at moderate microwave power P — 0.01 
[Eq. (8.2)] vs lo c /lj. Solid lines: the functions — d 3 [panels (a) 
and (c)] and h a [panels (b) and (d)]; (a) and (b) correspond 
to a-i/2 = 1, while (c) and (d) to Q1/2 = 10. Dashed lines: 
the linear-in-P asymptotes d s + d s and h a [Eqs. (5.15), 
(5.16), and (5.21)]'. 
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FIG. 8: Solid lines: same as in Fig. 7 for P = 1. Dashed lines: 
the asymptotes for d s and h a in the SIC regime [Eqs. (8.5) 
and (8.7)]. Dotted lines: the saturated inelastic contribution 
[Eq. (6.6)]. 
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FIG. 9: Same as in Fig. 8 for P = 10. 



FIG. 7: Solid lines: same as in Fig. 6 for P = 0.1. Dashed 
lines: the asymptote for h a in the SIC regime [Eq. (8.7)]. 
Dotted lines: the saturated inelastic contribution [Eq. (6.6)]. 



[Eq. (5.21)]. The displacement part d[ A ^ is small due to 
the large ratio Ti n /r q = 50, while the quadrupole con- 

tribution to the Hall part, h s [Eq. (5.20)], as well as 
the anisotropic diagonal part (5.17) are absent for cir- 
cular polarization. In the case of linear polarization, 

(c) 

h s would provide a Hall contribution of the same or- 
der of magnitude as ha . For comparison, the linear 
asymptotes di B ^ + di D ^ and hi°^ are shown in Fig. 6 by 
the dashed lines. The mismatch between the solid and 
dashed lines is negligible for small detuning from the res- 
onance, 



A = \uj/uj c - 2| < 1 , 



(8.4) 



owing to the sin 2 7rA term in Eq. (8.1). Away from the 
resonance, at A ~ 1/2, the parameter Q ln ~ 0.1 and the 
feedback effects (Sec. VI) lead to a noticeable correction 
to the linear behavior. 

The evolution of d s and h a with increasing P reveals 
the nonlinear effects studied in Sees. VI and VII. First of 



all, note that the part of <7 p h associated with Fid and Fi^ 
[Eqs. (7.25) and (7.26)] vanishes both at A = and at 
uj c /oj = 2/5 and 2/3, where A = 1/2. The photoresponsc 
at both ends of the interval is, therefore, fully governed 
by the contributions to <r p h coming from F2d, F^d, F2h 
and F3/j, which represent subleading terms for small P. 
In accord with the results of Sec. VII, Figs. 7(a) and (c) 
show that the photoresponse at A = 1/2 in d s is indeed 
discernible only at P > a -1 ; otherwise, it is masked by 
the much stronger inelastic contribution that develops at 
A<l/2. 

In the vicinity of the resonance, the diagonal photore- 
sistivity is dominated by the Fid term in Eq. (7.25), 
which at A 2 <C (Pa) -1 can be represented as 



d s ~ w cot w Fid 



2ttoj .P;* sin(27ra;/u; c ) 
to c 1 + sin 2 (ttuj J ui c ) 

P sin 

2uj c uj c 

TTUJ 2tTLU ( P* TTUJ 

- — sin sin — 

AUJ C UJ C \0JT q UJ C 

{uj + uj c ) 2 + 3w 2 



{uj + 2uj c ) 2 



(8.5) 
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where 



(uj + uj c ) 2 



p* = ln p* 



(8.6) 



see Eqs. (6.18), (6.20), and (7.29). The asymptote (8.5) 
is shown in panels (a) and (c) of Figs. 8 and 9 by the 
dashed line, while the dotted line describes the behavior 
of the first term of Eq. (8.5), which represents the satu- 
rated inelastic contribution [Eq. (6.6), see also Eq. (15) of 
Ref. 50]. For P ^> aT 1 the photoresponse tends to con- 
centrate at A <C 1, which can be seen from its evolution 
with increasing P in Figs. 6-9. 

Similarly, the Hall part at A 2 <C (aP) -1 can be ap- 
proximated by 



K - 



2u c P* sin 2 (7rw/w c ) 



U)T q (uj + 2bJ c ) 



+ 2tt 1 



P* Sm 2 (TTtd/uJ c ) 



(uj + 2cj c )r q 
nP* sin(27ro)/w c ) 1 - P* n sm 2 (Trw/u c ) 
(cj + 2uj c ) Tq l + P* n sin 2 (7ru>/w c ) ' 

(8.7) 

see Eqs. (6.10), (7.31), and (5.10). The asymptote (8.7) 
is shown in panels (b) and (d) of Figs. 7-9 by the dashed 
line. 

Let us now briefly recall the physics behind the cal- 
culated dependencies by the example of Fig. 9, where 
the microwave power P = 10 is large and the effective 
power Q [Eq. (8.1)] changes from Q — at u> c /u> =1/2 
to Q ~ 10 at uj c /uj — 2/5 and 2/3. As a result, in Fig. 9, 
the system passes with increasing A through all regimes 
discussed in Sees. V-VII. 

• A = 0: exactly on resonance, the photoresponse is 
zero. 

• A -C A max <C 1: the linear photoresponse, 
d s , h a oc PA (Sec. V). The inelastic contribution 
dominates the photoresponse in the diagonal part, 

4 = -^(2^) 2 — PA (8.8) 

9 Tq 

[linear-in-P part of the first term in Eq. (8.5)], 
while in the Hall part the photovoltaic mechanism 
produces 



h a = --(2irfa 1/2 PA 



(8.9) 



[linear-in-P part of the last term in Eq. (8.7)]. The 
maximum value of the inelastic term in d s as a func- 
tion of A is reached at 



A max = Y / 9r q /2 7 r 2 T in P 



• A <~ A max : the crossover to the SIC regime 
(Sec. VI). As detuning from the resonance A ap- 
proaches A max , the nonlinear corrections in the 
first term of Eq. (8.5) and in the last term of 
Eq. (8.7), driven by the feedback effects related to 
the strong oscillatory pattern in the isotropic part 
Poo (e) of the distribution function, become impor- 
tant. At A = A max , the inelastic contribution to 



d s reaches maximum, d s = ±4A 



\/P, and 



(8.10) 



crosses over into the A -1 decay at A > A max (il- 
lustrated by the dotted lines). In the Hall term 
h a , the strong oscillations of Poo change sign of the 
photovoltaic contribution at A = A max . 

• A max < A < {Pay 1 / 2 : the SIC regime, in which 
the photoresponse becomes independent of T; n and, 
hence, on temperature. This property of the high 
power photoresponse can be easily verified experi- 
mentally 

At A ^ A max , the two last terms in Eq. (8.5) be- 
come relevant, leading to a decrease of the diagonal 
photoresponse. The last term is actually irrelevant 
in Fig. 9(a) (which corresponds to a\j 2 = 1), so 
that [Eq. (7.29)]: 

<i s c— cot— (4-9Q)^|-8tt 2 PA. (8.11) 

The second term in Eq. (8.11) consists of the dis- 
placement contribution modified by the strong os- 
cillations in P o [in the SIC regime, its sign is in- 
verted compared to the linear case, Eq. (6.12)] and 
the two-photon correction to Poo itself [Eq. (5.3)]. 
This term is responsible for the strong narrowing of 
the oscillation in d s as a function of A and also for 
the nonmonotonic dependence on power at fixed A. 
At a » 1, the last, photovoltaic term in Eq. (8.5) 
becomes important, leading to additional narrow- 
ing of the peak in Fig. 9(c) and Fig. 8(c) compared 
to Fig. 9(a) and Fig. 8(a). 

The Hall part h s exhibits an additional oscillation 
at A <~ (Pa) -1 / 2 and decays for larger A due to the 
feedback and multiphoton effects. At A ^ A max , 
the first two terms in Eq. (8.7) (which, in contrast 
to the last term, are even in detuning lo — 2lo c from 
the resonance) become pronounced, as can be seen, 
for example, from the difference in amplitude of the 
second peak and the second dip in Fig. 9(d). 

• A > (Pa)- 1 / 2 : the ultrahigh power regime, m 
which all effects related to the Landau quantization 
decay to zero due to the feedback and multiphoton 
effects. 

Figures 6-9 demonstrate that, for large P and/or large 
a, the strongly nonlinear effects studied in Sec. VII es- 
sentially modify the behavior characteristic of the SIC 
regime 50 (dotted lines in Fig. 7-9). In particular, high- 
power measurements should reveal the significant nar- 



20 



rowing of the oscillation in p xx with increasing P, mani- 
fest in Figs. 8 and 9. It should be noted that, in addition 
to the magnetoresistivity traces dominated by the strong 
feature in the vicinity of A = 0, it is worth studying the Q 
and T-dependences of the photoresponse at several fixed 
ratios w/w c ; in particular, the contributions ¥24 and 
which manifest themselves at A <~ 1/2. 

At P ~ or 1 the magnitude of the photoresponse 
in p xy , Fig. 7(d), becomes comparable to the effect in 
p xx , Fig. 7(c), which seems to be in accord with the 
experiment. 13 However, the experiments that revealed 
the Hall oscillations with h a ~ — d s were performed at 
ujT q <~ 1, i.e., at a ~ 1, where for the level of power 
used in the experiment the predicted oscillations in h a 
should be more than an order of magnitude smaller than 
those in d s . Our theory is not capable of explaining the 
observed effect in p xy (probably related to microwave- 
induced u) /w c -dependent changes in the electron concen- 
tration, see discussion in Sec. VIII of Ref. 50). 



IX. RANGE OF APPLICABILITY OF THE 
THEORY 

The above calculation the OPC is based on the as- 
sumptions that (i) disorder is smooth, Eq. (2.2); (ii) tem- 
perature is high, Eq. (3.3); (iii) Landau levels strongly 
overlap with each other, Eq. (4.1). Altogether, these 
conditions made the analytical treatment of the problem 
possible to all orders in the microwave power. Now let us 
briefly discuss the range of applicability of the above ap- 
proximations, especially in the limit of strongly nonlinear 
photoresponse. In particular, it is important to discuss 
the role of heating. 



Here T e _ 1 ph (To) oc Xj 3 is the inelastic relaxation rate due 
to scattering on acoustic phonons (the latter represent a 
thermal bath at temperature T ) and s± describes po- 
larization of microwaves [Eq. (2.20)]. At a typical value 
of To ~ 1 K, the electron-phonon scattering rate exceeds 
t^ 1 by more than an order of magnitude, so that the 
parameter 77 in Eq. (9.2) can vary between 1 and 1000 
under realistic experimental conditions. At first glance, 
the large values of r\ imply a strong heating, T e » T . 
However, this is not actually the case in view of the 
strong temperature dependence of T _ p h(T). While at 
77 < 1 the temperature grows linearly with 77 (specifically, 
T e — To ~ 77), in the limit of strong heating the growth is 
much weaker, T e /T ~ 77 1 / 5 3> 1. The weak dependence 
on the microwave power at large 77 means that, in prac- 
tice, the heating can be substantial but never strong, 
i.e., T e /To ~ 2 3. At the same time, with lowering 
bath temperature well below 1 K, while keeping the mi- 
crowave power constant, the heating effect may become 
significantly more pronounced (since 77 grows with lower- 
ing temperature as T Q ~ ), in effect restoring the high-T 
limit we deal with in the paper. 

In the dc response, the strongest manifestation of 
the heating effect is the exponential suppression of the 
Shubnikov-de Haas oscillations (provided the latter are 
visible in the absence of microwaves at T e = To). The 
OPC is influenced much more weakly, since T e enters 
the photoconductivity only indirectly, through the T- 
dependence of the electron-electron inelastic scattering 
rate, t^ 1 oc T e 2 . In particular, in the linear regime, 
P <C T q /ri n , taking the heating effects into account leads 
to a sublinear dependence of the diagonal photoresponse 
on the microwave power, namely d s oc PT~ 2 (P), while 
at high power, P » Tq/im, the maxima and minima of 
d s scale as s/P/T e {P). 



A. Heating effects 



B. Stratification of the distribution function 



The high-T approximation led us to Eq. (3.3) for the 
oscillatory distribution function (see Appendix A). The 
Fermi-Dirac part of the distribution is characterized by 
the effective temperature T e 3> u) c . The oscillatory part 
appears due to the Landau quantization, which, in the 
limit of strongly overlapping Landau levels, is weak: the 
amplitude of the oscillations of the DOS is proportional 
to the small parameter 5 [Eq. (4.1)]. The heating of elec- 
trons, which leads to a growth of T e with increasing mi- 
crowave power, can therefore be treated separately within 
the framework of the classical Boltzmann equation with 
a constant DOS. Such a calculation leads to the following 
equation for T e : 41 



(T e /T ) 4 (T e /T - 1) = 77 , 
_ Tc-ph(To) / ev F E u 

V ~ Ttr V 2T 



E 



U) ± U) c 



(9.1) 



(9.2) 



Apart from the condition on the temperature T e 3> 
1/ts, the approximation (3.3) is valid only if the os- 
cillatory correction to the distribution function remains 
much smaller than unity, so that the general requirement 
< /(e) < 1 is met. Despite the DOS being only weakly 
modulated (5 <C 1), the above condition is violated at suf- 
ficiently high power for the values of the ratio u/lo c corre- 
sponding to the minima and maxima of d s [see Figs. 9(a) 
and (c)] . Near these points, the effect is governed by the 
inelastic mechanism, Eq. (8.11), and the oscillatory part 
of Eq. (3.3) can be estimated as 



^(eR/ T 



Sd s 



sin^s) 



< 



t B T e A E 



(9.3) 



where we used Eqs. (4.3), (4.7), (6.4), (6.6), and (8.11). 

At large power, the maximum value of d s grows as 
y/P/T e (P) [see Eq. (8.8) at A = A max ]. It follows that at 
P <~ Pi, = (u) c S/T e ) 2 Ti n /Tq the oscillations of the isotropic 
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part of the distribution near the maxima and minima of 
d s become of order unity. At small detuning from the cy- 
clotron resonance harmonics, A ~ A max , and at P > P+, 
the splitting of /(e) into the smooth and oscillatory parts 
(3.3) is no longer possible. The distribution function in 
this regime can be calculated numerically according to 
Eq. (11) of Ref. 50 (for details of the numerical procedure, 
see Ref. 68). An illustrative example of the stratification 
of the distribution function at P ~ P* is shown in Fig. 10. 
Both the heating and the stratification of the distribution 
function suppress the growth and the narrowing of the 
peak in d s at very high power. Beyond the close vicinity 
of maxima and minima, at |A| = \lo/lo c — N\ 3> 5/tBT e 
(where N is integer), the oscillatory correction to the dis- 
tribution function remains small at any power level and 
the high-T approximation works well, so that our results 
remain intact whatever the microwave power. 



to 

0.5- 








e 



FIG. 10: Isotropic part of the distribution function at three 
levels of the microwave power: P = 0, P <C P*, and P ~ P*. 



C. Smooth vs short— range disorder 

In this paper we assumed that the disorder potential 
is created solely by charged impurities separated from 
the plane of the 2DEG by a spacer of width £/2 S> kp 1 
[Eq. (2.2)]. Such a minimal model of disorder allowed us 
to account for the experimentally relevant small-angle 
scattering condition r q 3> T tr and to consider in an un- 
ambiguous way the region of magnetic fields t^. 1 <C lo c <C 
t~ . On the other hand, a "two-component model" of 
disorder, 76 including strong scatterers on the background 
of the smooth potential (2.2), is known to provide a bet- 
ter description of real ultrahigh-mobility structures. 

To our current understanding, the inclusion of the 
strong scatterers into the theory should not lead to any 
qualitative change of the results presented here [provided 
the conditions (2.1) remain satisfied]. In particular, the 
derivation of the linear-order inelastic contribution (5.15) 
in Ref. 34 does not require any assumptions as to the type 
of disorder. 



By contrast, the strong scatterers were shown 3,66 to 
play a crucial role in the nonlinear dc response (in the 
absence of microwaves) , where strong oscillations of the 
magnetoresistivity ("HERO") were observed, 3,20 ' 21,28 ' 29 
governed by the ratio of the Hall electric field to the 
magnetic field. Remarkably, a strong interplay be- 
tween two types of oscillations in the 2DEG driven by 
both microwave and strong dc fields was recently ob- 
served experimentally. 30 Theoretical description of this 
experimental situation necessitates the inclusion of the 
backscattering off strong impurities, which we relegate 
to future work. It was argued in Ref. 57 that, in the 
regime of strong dc field and in the presence of short 
range disorder, effects of the microwave field on screen- 
ing of impurities should also be taken into account. 



D. Stronger magnetic fields, cu c T q > 1 

The calculations in Sees. V-VII were performed for the 
case of strongly overlapping Landau levels, cj c T q <1. In 
this limit, (i) microwave-induced corrections to the DOS 
can be neglected; (ii) there are substantial effects related 
to excitation of higher angular and temporal harmonics of 
the distribution function; in particular, the photovoltaic 
and quadrupole contributions to the OPC. At stronger 
magnetic fields, the modulation of the DOS due to the 
Landau quantization becomes more pronounced: in the 
SCBA approximation, at uj c T q ~ 1.8 the Landau levels 
get separated from each other. 68 For the case of separated 
Landau levels, the effect on higher harmonics of the dis- 
tribution function becomes strongly suppressed and the 
diagonal part of the photoresponse is determined by the 
interplay of the inelastic and displacement mechanisms. 

In addition to the nonlinear effects considered in this 
paper, at uj c T q > 1 the DOS itself is modified by mi- 
crowave radiation. In particular, microwave-induced 
sidebands appear on both sides of the Landau levels, 
leading to the appearance of an additional structure in 
the OPC pattern near the fractional harmonics of the 
cyclotron resonance, nui = muj c . 7S A similar structure 
in photoresponse also arises as a multiphoton correction 
to the inelastic contribution to the OPC. Both struc- 
tures have precisely the same shape (as a function of to); 
however, the multiphoton mechanism dominates in the 
case of strongly separated Landau levels. 78 Such a "frac- 
tional oscillatory pattern" (FOP) in the nonlinear pho- 
toresponse was observed in recent experiments 10 ' 26,27 and 
explained in terms of either of multiphoton corrections to 
the displacement contribution 10 ' 26,64 or, alternatively, in 
terms of single-photon resonant corrections to the inelas- 
tic contribution which are specific to the crossover region 
uj c T q ~ l. 27 > 68 All four contributions to the FOP will be 
compared and analyzed elsewhere. 78 
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X. CONCLUSIONS 

Summarizing, we have developed the systematic ap- 
proach to the microwave-induced oscillations in the mag- 
netoresistivity of a 2DEG. This approach has enabled us 
to classify contributions to the photoresistivity according 
to the combined action of the microwave and dc fields 
on the temporal and angular harmonics of the distribu- 
tion function (Sec. III). We have studied the interplay of 
the resulting mechanisms of photoresponse at high mi- 
crowave power. In the limit of strongly overlapping Lan- 
dau levels (Sec. IV), the dc photoconductivity has been 
calculated to all orders in the microwave power (Sec. V- 
VII). 

To linear order in the power, two novel mechanisms of 
the oscillations (quadrupole and photovoltaic) have been 
identified, in addition to the known ones (displacement 32 
and inelastic 34 ). The quadrupole and photovoltaic mech- 
anisms have been shown to be the only ones leading 
to oscillations in the Hall part of the photoconductiv- 
ity tensor. Of particular interest is the result that the 
quadrupole contribution violates Onsager symmetry. In 
the diagonal part, the inelastic contribution dominates at 
moderate microwave power, while at elevated power the 
contributions of other mechanisms become important. 

In Sees. VI and VII we have considered the strongly 
nonlinear photoresponse at high microwave power. We 
have shown that a competition between various nonlin- 
ear effects (the feedback effects, the excitation of high 
angular and temporal harmonics of the distribution func- 
tion, and the multiphoton effects) drives the system 
through four different nonlinear regimes with increasing 
microwave power. 

Most dramatic changes in the photoresponse are due 
to the feedback effects. In the SIC regime, Sec. VI, the 
feedback from the microwave- induced oscillations of the 
isotropic part of the distribution, F o, leads to the satu- 
ration of the inelastic contribution, 40,50 and to the strong 
interplay of the inelastic effect and all other contributions 
to the OPC. In particular, the strong oscillations of i*oo 
change sign of the most relevant parts of the displacement 
and photovoltaic contributions. Remarkably, in the SIC 
regime the photoresponse becomes independent of the 
inelastic scattering rate and, therefore, of temperature. 

At higher power, Sec. VII, the feedback suppresses the 
effects on higher temporal and angular harmonics of the 
distribution function. At still higher power, the mul- 
tiphoton excitation becomes pronounced and starts to 
compete with the feedback effects. At ultrahigh power, 
the feedback and multiphoton effects destroy all quantum 
contributions, restoring the classical Drude conductivity. 

Our theory predicts nonmonotonic behavior of the di- 
agonal photoresponse as a function of the microwave 
power at a fixed ratio oj /uo c . As illustrated in Sec. VIII, in 
magnetoresistivity measurements such a nonmonotonic 
dependence should result in a significant narrowing of the 
oscillatory structure around the integer values of lo/lu c 
at sufficiently large power, see Eq. (8.11). To a large ex- 



tent, the narrowing is due to the displacement contribu- 
tion. At the same time, the amplitude of the oscillations 
is controlled by the inelastic mechanism at any power 
level. We suggest to experimentally measure the power 
and temperature dependences of the photoresponse at 
fixed u)/u) c not too close to the cyclotron resonance har- 
monics. Such experiments will make it possible to ob- 
serve the rich behavior of the photovoltaic contributions. 
Also, our theory predicts oscillations in p xy , which, in 
the limit of strongly overlapping Landau levels, can be 
comparable in amplitude with those in p xx . 
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APPENDIX A: QUANTUM BOLTZMANN 
EQUATION IN THE HIGH-TEMPERATURE 
LIMIT 

In this Appendix we derive the high-T version of the 
quantum Boltzmann equation starting from its general 
form, Eqs. (2.11), (2.13). We show that the solution of 
the QBE in the limit T 3> lo c has the form 



/ 21 = /(t_,i)=0(i_)+ J2 *n(V,t)i«(tL n) ) 

Tl— — CO 

(A.l) 

and obtain the equations for the amplitudes <&„(<£, t). 
Here t_ = t 2 - h, t = (t 2 + h)/2, and = 
(27r) _1 Jde /t(e) exp(— iet-) is the Fermi distribution 
function in the time representation, 

»rexp(-fe F i-) (n) 
y(t ~ } 2sinh^T(t-+zO) ' " *" B ■ 

In the energy representation, the solution (A.l) ac- 
quires the form (3.3), 

/(e,¥>, i) - jdt-e ist -f 21 

= f T (e)+F(e,t,<p)d £ Me) , (A.2) 

oo 

T{e,Lp,t) = ^2 ®n{<P,i) exp(inet B ) ■ (A. 3) 

n— — oo 

In the high-T limit, the function /(e) shows fast oscil- 
lations on top of the smooth thermal distribution /t(e)- 
The function / 2 i, Eq. (A.l), is composed of equidistant 
sharp peaks of width f/T, with the distance between 
them being equal to ts- 
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The DOS (3.1) for high Landau levels is a periodic In terms of the coefficients g m , the QBE is rewritten as 
function of energy and the functions g R — — (g A )\ en- 
tering Eq. (2.13), obey + Uc d v )f 31 - St in {/} 3 i = St im {/} 31 (A.5) 



m=0 



(A.4) with 
I 



Stim{/}31 



E 

m>0 

E 

m>0 



i h + ts\ ~ (t 2 + h 

^3l9m 7, Hi - /21/V-32ffr 



, , ,. ,h + h\ ~ „ (t 2 + ti 

^31 J 329m = - /32Ml5 r 



t 2 =t 3 -mt E 



t2=*l —mt B 



(A.6) 



Now we substitute the solution in the form (A.l) into the kinetic equation. The unperturbed part G(t~) gives zero 
on the l.h.s. of Eq. (A.5). The impurity collision integral acting on G(t~), produces 



St im {£} 3 i = Q ^ ~ ~ ^9m{t 3 - mi s /2) 



m>0 



t3—t2=mtB 



ti —t 2 =mt B 



(A.7) 



where we took into account the fact that the operator /C, in the time representation, commutes with an arbitrary 
function independent of ip, sec Eq. (2.25). 



Let us consider the first sum in the above expression. At 
TtB 3> 1, the (5-function-likc G(t 2 — ii) puts t 2 within the 
interval \t 2 -h\<l/T. Strictly at the point t 2 -h = 0, 
the function Q(t 2 — t\) is infinite, while the expression 
in the square brackets is zero. Recalling that the kernel 
changes smoothly on the scale of \t 2 — t\\ ~ 1/T <C ts, 
one can replace the difference in the brackets by the 
derivative d\JCzi- The result is a series of peaks of iden- 
tical shape, 



An>o = 
An<o = 
A N =o = -2i 



[(d t _ -d t /2)t 21 
(d t _ + d t /2)t 21 

d t _fc 2 i 



t-=Nt B 



9N(t) 



t_=0 



t_=ATt B 

.go = -i d t _JC 



9- N (t) 



t_=o 



(A.8) 



where the square brackets mean that the time derivatives 
act on the kernel only and do not act on the functions 
gN(t)- In the last expression, we took into account that 
go, representing the DOS in the absence of magnetic field, 
is not affected by the external fields and by disorder: the 
average DOS is a conserved quantity, go = 1/2- 

One can see that Sti m {<7}, given by Eq. (A.8), has ex- 
actly the same form as the oscillatory correction to the 
distribution function, Eq. (A.l). It follows that the solu- 



tion at Tts ^ 1 indeed has the form (A.l). Substituting 
the perturbed part of the distribution function (A.l) into 
the kinetic equation, Eqs. (A.5), and (A.6), we finally ar- 
rive at the equations for the amplitudes $„, 

(d t + Lu c d v )<P N + Tr 1 ^) = A N 

+ |^iv(*)$ 'N -m(t m )9m(tm-N) 
m>0 

— &N-m(t m )lC m (t m -N)gm(tm~N) 
+JC N (t ) $ N+m {t m )9m {t m + N ) 

— <&N+m(tm)fc—m(t m +N)gm(tm+N) , (A. 9) 

where IC n (t m ) = IC(t- — nts, t = t m ) and t m — 
t — mtB /2. The inelastic relaxation, which controls the 
magnitude of the oscillations in the isotropic part of the 
distribution function, can be described in the relaxation 
time approximation, with T ln being the effective electron- 
electron inelastic scattering time. 50 



APPENDIX B: QUANTUM BOLTZMANN 
EQUATION: OVERLAPPING LANDAU LEVELS 

In this appendix we consider the high-T version of 
the QBE, Eq. (A. 9), in the limit of weak magnetic field, 
which corresponds to the exponentially small modula- 
tion of the DOS. In this case, only first two terms in the 
expansion (A.4) should be taken into account, namely, 
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go = 1/2, and gi = —S, where 5 = exp(—t B /2T q ) <C 1. 
Other terms, g n = 0(S n ), are exponentially smaller and 
can be neglected. We will also use the smallness of the 
odd part of the kernel, JCj, with respect to JC±, con- 
trolled by the parameter y/ Tq/rtr- At zero order in JCj, 
the distribution function is even in <j) and the current is 
zero. Thus, JCj -C JC± should be taken into account once, 
while the contributions of the second and higher orders 
in JCj should be neglected. 

At leading order in 5, the kinetic equation (A. 9) gives 



(d t + <jj c d v )$ = -i d t _JC 



= —i d t _iCj 



(B.l) 

where the operator JC at the last position in any expres- 
sion should be understood as acting on unity. The in- 
elastic term on the l.h.s. is omitted as the resulting $ 
contains odd angular harmonics only. Indeed, according 
to Eq. (2.21), X(t- = 0) = 0. It follows that only odd 
part of the kernel, JCj, produces a non-zero contribution 
to the r.h.s. of Eq. (B.l) [see Eqs. (2.25)-(2.28)]. Also, 
on the r.h.s. of Eq. (B.l) we omitted 



[£o(t)$o(t) - *o(*)£o(*)]/2 



--1^2 



(B.2) 



originating from the m = term in Eq. (A. 9). (Here 
we used again the fact that X(t- = 0) = 0, which gives 



/Co 



1 + JC, 



). In the case of classically strong 



magnetic held, w c r tr 3> 1, the term (B.2) can be safely 
neglected. Note, however, that the dc conductivity is 
infinite in the absence of this term in the opposite limit, 
5 = 0. 

At order O(S), the kinetic equation (A. 9) produces an 
oscillatory correction with the amplitude $i(tp,t) obey- 
ing 

(d t +w c v )$i(t) + r i ; 1 ($ 1 ) = iS(d tB - d t /2)t(t B ,t) 

+ [JC 1 (t) - £,] *i(t) - 6 [£i(i)* (ti) - *o(ti)£i(*)] • 

(B.3) 

In classically strong magnetic held, the last term should 
be neglected, since (i) $o, according to Eq. (B.l), is gen- 
erated by JCj oc Tq/rtr -C 1; (ii) the leading part of the 
kernel, JC±_ oc (r q /T tr ) , commutes with ip. Thus, only the 
odd part of the kernel, JCj, produces a nonzero term in the 
curly brackets, and the result, of second order in JCj, is 
proportional to r q /r tr and should be neglected. With the 
same accuracy, C(^/r q /rt r ), the expression in the square 

brackets should be replaced by [JCj(t) + JC±(t)] t _=t B ■ 

In the high-T limit, the leading contribution to the 
current is of order 6 2 [Sec. IV] . Accordingly, we represent 
the nonoscillatory part of the solution (A. 3) as 



with (j>Q D ^ obeying Eq. (B.l) and (j)^ generated by m = 1 
term in the sum (A. 9), 

25(d t +io c d v )^ ) (t) 

= -^ (t)$_ 1 (t 1 ) + *_i(ti)^i(ti) 

-£o(t)*i(ti) + *i(*i)^-i(*i) ■ (B.5) 

Now we recall that the distribution function (A. 3) is real, 
= JC n = lC- n , and take into account that Im/C = 
—iJCj. After that, Eq. (B.5) can be rewritten as 

(dt+UJcd^W = 

Im[0i(fi)]z/Cj(i s ,ti) - Re[<j) 1 (t 1 )]JC±(t B ,t 1 ) , 

(B.6) 

where the part JC V of the kernel is neglected, t\ = tg/2, 
and we switched to the more convenient notation of 
Sec. IV by putting <j>i = -$i/<5. 

Equation (B.3), being rewritten separately for the real 
and imaginary parts, reads 

= (d t /2 ~ d tB )JC ± (t B ,t) + JC ± (t B ,t)lmMt) , 

(B.7) 

(d t + cu c d v )RcMt) - (dt/2 - d tB )itj(t B ,t) 
+t^(t B ,t)Re<j) l (t) + itj^B^lm^t) . (B.8) 

Equation (B.8) shows that Re0i is of hrst order in JCj, 
hence an odd function of ip. Accordingly, we omitted the 
inelastic collision term in Eq.(B.8), as well as the term 
iJCj(t B ,ti)Recj)i(t) on the r.h.s. of Eq.(B.7) which is of 
second order in JCj and thus small in parameter T q /r tr . 

The similarity of the r.h.s. of Eq. (B.8) to that of 
Eq. (B.6) is not accidental. In fact, the dc current is 
expressed in terms of the combined quantity, <fij(t) = 

</>o 2) (£ + W2) + FW>i(t), see Eq. (4.8). Both Eqs. (B.8) 
and (B.6) contain the term JC±(t B , t)Ke(f)i(t) = 
Re(pi(t)JC±(t B ,t) [unlike JCj, the operator JC± commutes 
with ip, see Eqs. (2.26), and (2.27)]. This term gener- 

(2) 

ates contributions to <p and Re0i(f), as illustrated by 
diagrams (E) and (F) in Fig. 1. However, in the dc cur- 
rent (4.8) these contributions cancel each other, as the 
r.h.s. of the equation for tfij contains JCj only, whereas 
4>±_ (t) — Im0i (t) is fully governed by JC ^ . The hnal equa- 



$0 = ^+2^ 



>o 2) 



(B.4) 



tions for 
Sec. IV. 



and <j>j, Eqs. (4.9)-(4.11), are given in 
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